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A cavity perturbation method for wide temperature range 
dielectric measurements is described. It is used to measure 
the complex dielectric constant of aqueous-DMSO solutions in 
the temperature range -60°C to +100°C, for concentrations 
from 0 to 100%. The results indicate that the relaxation 
frequency drops from 18 Ghz to 3 Ghz as the concentration is 
increased from 0 to 70%. A further concentration increase 
Causes the relaxation frequency to increase to 9 Ghz. At 
temperatures below freezing, the value of €" is increased by 
the addition of DMSO, and its rate of change with 
temperature is reduced. The cryobiological implications of 


this result are examined. 


A numerical model of plane wave microwave heating of 
Canine kidneys is developed. Details of 1.5 mm and larger 
are represented in the model. The model shows that microwave 
energy absorbtion is enhanced by DMSO perfusion, and that 
adjustment of DMSO concentration can provide control of 
heating uniformity. Uniform perfusion is found to be 
essential if heating is to be uniform. Steep gradients in 
the electric field are found to exist at the interfaces to 
the renal sinus. The compositions of the solutions residing 
in the renal peivis and vascular structure of the kidney are 


found to have considerable effect on the heating uniformity. 
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Impedance of free space 
Temperature 

Initial temperature 

Final temperature 
Wavelength in tissue type on 
Free space wavelength 
Maximum wavelength in scatterer 
Permeability 

Permeability of free space 
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(1.1) Microwave Heating. 


The expansion of the biological sciences has led to a 
need for techniques for the preservation of biological 
material. Low temperature storage, as developed by the 
science of cryobiology, has provided one solution to this 
problem. There are three phases involved in the process of 
low temperature storage. These are: (a) the cooling of the 
tissue, (b) its storage at a low temperature and (c) its 
recovery from the frozen state. The present work is 
primarily concerned with the use of microwaves in the third 
phase, though certain aspects of the first two phases are 
important. Microwave heating is dependent on the dielectric 
constant and so factors which influence this constant must 
be considered. In particular, account must be taken of the 
type of cryoprotective agent, the degree of perfusion with 


the agent, and the temperature range. 


The use of microwaves for the recovery of frozen organs 
was first reported by Andjus (1) in 1953. Subsequent 
investigations, including comparisons with conduction 
heating, have been made by Kenney (2), Holst (3), Hamilton 
(4), Ketterer (5,6), Burns (7) and others. The major 
difficulties encountered are achieving and controlling the 
heating rate, attaining reasonable heating uniformity, and 
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stopping the process at some desired temperature. Research 
aimed at resolving these problems has lagged considerably 
behind the other areas of cryobiological investigation, but 
the increaSing availability of microwave apparatus is aiding 
in these studies. Systems designed to solve one or more of 
these problems have been described by Rzepecka (8), Rajotte 
(9), Walker (10) and others, although most of these systems 


are primarily useful for small samples. 
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The two types of heating employed in cryobiology are 
microwave heating and conventional heating, which employs 
heat transfer by conductive or convective techniques. In the 
latter case, it is usually necessary to immerse the tissue 
in a liguid (often water) to prevent surface evaporation. In 
microwave heating, however, the surrounding air temperature 
is not affected, and so liquid immersion is not necessary. 
If it is reguired, liquids with suitably low dielectric 


loss, such as fluorocarbon, can be used (11). 


For conventional heating, the delicate nature of 
biological materials limits the temperature of the 
surrounding liguid to less than 509°C. This effectively 
limits the heating rate. As the size of the heated sample 
increases, this limitation becomes more severe. in addition, 


heat transfer requires that there be a temperature gradient 
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from the sample surface to its interior, and so acceptably 
uniform heating is not possible except for very small 
Samples. These two factors represent a physical limitation 


of conventional heating. 


In the case of microwave heating, energy is dissipated 
throughout the sample. Heating uniformity is a function of 
the disposition of the dielectric constants in the sample 
and the design of the heating apparatus, factors over which 
there is some degree of control. In addition, heating rate 
is limited by the available microwave power. Therefore, 
uniformity and heating rate are governed by technological 


factors, and so the potential for improvement does exist. 


There are three main problems with microwave heating. 
These are: (a) The amount of dissipated energy is governed 
by temperature dependent dielectric properties. (b) 
Temperature measurement during heating is difficult because 
of electromagnetic field perturbation by the measuring 
probe. (c) The coupling structure must often be specifically 
designed for a given sample if optimum uniformity and energy 


transfer are to be achieved. 


An investigation into the temperature dependence of the 
dielectric constants of tissue and cryoprotective agents 
will aid in the solution of the first probiem, and this is 


dealt with in Chapter 4. A viable first order solution to 
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the second problem has been presented elsewhere (10). The 
humerical model of microwave heating dealt with in Chapters 
5 and 6 represents a first step towards the solution of the 


third problen. 


(1.3) Cryoprotective Agents and their Uses. 


When biological tissue is cooled through the 
temperature range from -60°C to +109C, ice-formation causes 
Changes in intercellular fluid concentrations with resulting 
Changes in osmotic pressure. These effects cause potentially 
damaging stresses to be established within the cellular 
structure of the tissue, and for this reason, the present 
microwave heating study is limited to this crucial 


temperature range. 


The cryoprotective agents used to protect tissue from 
damage during freezing are basic to the science of 
cryopreservation. There are two major types of 
cryoprotective agent commonly employed. These are: (a) some 
relatively low molecular weight materials (70-400D) such as 
dimethyl sulfoxide (DMSO) and glycerol, which act by 
permeating the intercellular spaces of the tissue, and _ (b) 
some high molecular weight materials such as polyvinyl- 
pyrrolidone (PVP) and dextran which are extracellular and 
are believed to lessen the rise in extracellular ice- 


concentration as ice forms, thereby protecting the cell 
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surface rather than the cell interior (12). The two types of 


agent are sometimes used in conjunction with one another. 


DMSO is the liquid which is primarily investigated in 
this work since it is both an effective cryoprotective 
agent, and, being a highly polar liquid, is of some interest 
in the area of dielectric theory. Ashwood-Smith (13) has 
reviewed many of the important properties of DMSO as related 
to cryobiology, and Szmant (14) has summarised some of the 


physical properties of the liquid (see Appendix C). 


The compound DMSO was first prepared in 1866, but only 
came into general use after World War II when its 
radioprotective and cryoprotective properties were 
discovered. Its ability to penetrate cell membranes aroused 
considerable interest in biological and medical circles, but 
reports of disturbances to the lens of the eye halted 
clinical tests. These tests have again begun, but the full 
potential of the compound has yet to be realised, though it 
has certainly proved important in cryobiology. The uses and 
properties of the liquid are discussed in considerable 
detail in the literature (15). Values of some dielectric 
properties are included there, while other measurements have 
been reported by Bourgoin (16) and Doucet (17), though these 
were not at the freguency and concentrations employed in the 


present work. 
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(1.4) Lhe Role of Dielectric and Thermal Properties. 


A considerabie amount of dielectric data for biological 
materials has appeared in the literature during the past 
three decades. Among the major contributors are England (18) 
and Schwann (19), with extensive compilations by Tinga (20) 
and Presman (21). Knowledge of these data are important for 
theoretical investigation in several areas. These are: (a) 
For analysis of electrocardiography and other related 
electrical measurements of body function. (b) For the 
detemination of the electrical properties of tissue at high 
frequencies. In microwave diathermy, this is a pre-requisite 
to understanding how high frequency electrical currents 
penetrate the body and produce heat. (c) For the analysis of 
the deleterious effects of electromagnetic radiation on 
humans. (d) For analysis of microwave cooking of foods and 
other industrial processes. (e) For theoretical studies 
involving the formulation of equations describing material 


behaviour over a wide range of frequencies and temperatures. 


In each case, the conditions under which the data is 
obtained depend on the intended application. For example, 
data for frozen foods as well as foods at room temperature 
are important in studying microwave thawing and tempering, 
whereas bulk properties of human body tissue at normal body 
temperature (37°C) are of interest in the study of the 


biological effects of microwave radiation. The theoretical 
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study of microwave heating in cryobiology is still in its 
infancy and so dielectric properties measured under the 
conditions pertaining to cryopreserved tissue are only just 
beginning to appear in the literature (22). These conditions 
include tissue which can be in a wide range of temperatures 
and is often perfused with a cryoprotective agent. A wide 
variety of cryoprotective agents can be used and it is 
difficult to estimate the degree of perfusion for any given 


CaSee. 


There is some evidence to suggest that the dielectric 
properties of tissue are governed by the properties of the 
interceilular fluids for frequencies in the microwave region 
(19). Further, there is an interest in the dielectric 
properties of aqueous solutions of the cryoprotective agent, 
DMSO, for a wide range of temperatures. On these grounds, it 
was decided to make the dielectric property measurements on 
DMSO for the conditions mentioned, and then, in conjunction 
with literature data for appropriate frozen foods, to obtain 
estimates for the dielectric properties of perfused tissues. 
These properties are required for the numerical microwave 
heating model. Some data for perfused kidney tissue is 
available in the literature (22), and this is used to aid in 


making the estimates. 


Although they are not used in this work, thermal 


properties are also important and a detailed examination of 
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microwave heating should consider the specific and latent 
heats and the conductivity of the material. Specific heat 
and conductivity are also temperature dependent properties, 
but there is little data regarding their behaviour for the 
materials of interest in cryobiology. Tressler (23), Awberry 
(24), Dickerson (25), Ohlsson (26) and Hill (27) have made 
measurements on the specific and latent heats and 
conductivities of various meats and other foodstuffs, and 
some conciusions regarding the behaviour of frozen organs 


can be drawn from these measurements. 


It is recognised that complex materials like biological 
tissues do not freeze or thaw at a single well-defined 
temperature, but rather over a range of temperatures. 
Therefore, to make heat transfer calculations more 
manageable, Dickerson (25) has suggested the use of the 
enthalpy of the materials of interest. The existence of 
tables of enthalpy for various bio-materials under various 
conditions would be a great asset to microwave heating in 


cryobiology. 


Dielectric and thermal properties are thus a pre- 
requisite for optimum design of microwave heating apparatus. 
In addition, reasonably accurate data regarding these 
properties are especially important when numerical modelling 


technigues are employed. 
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(1.5) The Importance of Numerical Modelling. 

Numerical modelling techniques have the advantage of 
being easy to apply when many different variables are 
involved. In the case of heating studies, the effects of 
changes in dielectric and thermal properties can quickly be 
assessed, and possible solutions to problems of heating 
uniformity can be assessed In many Situations, the 
measurements required to characterise system behaviour are 
difficult to perform, and the use of numerical modelling is 


the only remaining recourse. 


FOr instance, there are considerable difficulties 
involved in making temperature measurements on microwave- 
heated materials. These difficulties are compounded if 
temperature distributions are required, as is the case if 
heating uniformity is of importance. Furthermore, the 
existence of numerous variables such as the composition of 
the cryoprotective agent, the heating rate adopted and the 
possibility of a combination of thermal and microwave 
heating mean that the direct experimental approach, if at 


all possible, would require many replications. 


The use of a numerical modelling techniques in this 
Situation would allow the interdependence of the different 
variables to be examined, and could lead to new techniques 


for controlling uniformity. One possibility would be the 
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introduction of liquids with carefully chosen dielectric 
properties to specific areas of the sample being heated. A 
Similar approach using metal balls rather than liquids has 
been employed in the literature (7). Another possibility is 
the simultaneous use of microwave and conduction heating, or 
of multifrequency heating (7). The evaluation of these and 
other approaches could certainly be simplified by the use of 


modelling techniques. 


In the past, theoretical microwave heating studies have 
been performed on geometrically simple shapes such as 
cylindrical *roasts' (28) and infinite slab ‘hamburgers'* 
(29). If inhomogeneous dielectric properties were considered 
at all, they were usually arranged in a geometrically simple 
fashion. For these reasons, the bulk properties of the 
materials of interest were usually sufficient. In the 
present work, a numerical model which is able to account for 
relatively small (1.5mm) structural inhomogeneities is 
employed, and so the properties of the individual tissue 
types with and without perfusion over a range of 
temperatures are reguired. More detailed models 
investigating other microwave effects have appeared in the 
literature (30,31,32), but body temperature properties 
(+109°C) of normal tissue were usually required in these 


Cases. 


The model adopted in the present work considers plane 
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wave illumination of a canine kidney. Though some situations 
justify the use of a plane wave, the model is not restricted 
to this type of iilumination. It is possible to simulate 
different types of incident wavefront, and even to model 
heating within a cavity. The main limitations are memory 


capacity and processing time. 


Chapters 2 and 3 deal with the development of a system 
to measure dielectric properties over a wide range of 
temperatures. Some of the problems encountered with the 


measurements are discussed in Chapter 3. 


In Chapter 4, the results of the dielectric 
measurements are presented, and some of the theoretical and 
biological implications of these results are discussed. In 
particular, the behaviour of the water-DMSO system in 
relation to its phase diagram is examined. In addition, 
Changes in the dielectric properties as a function of 
temperature are discussed in relation to their effect on 


microwave heating. 


Chapter 5 deals with the creation of a numerical model 
designed to study microwave heating effects. Application of 
this model to the canine kidney is described in Chapter 6, 


where the results obtained with the model are also 
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discussed. 


The two sections of Chapter 7 deal respectively with 
Suggested further research on the measurement of dielectric 


and thermal properties, and the use of the numerical model. 
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There are a great many different methods of making 
dielectric measurements, and it is not the present purpose 
to exhaustively review these techniques. The subject of 
dielectric measurements in general is adequately treated in 
the literature, for example, Bussey (33) and Altschuler 


(34). 


Broadly speaking, there are three general classes of 
dielectric measurement technigue for use at microwave 
frequencies. These are (a) Standing wave methods, (b) 
Perturbation methods and (c) Time domain techniques. For a 
given application, individual methods are judged according 
to the following criteria: 

(1) Accuracy and precision of the method. 

(2) Sample size and configuration. 

(3) Frequency range where the method is applicable. 
(4) Available and required apparatus. 


(5) Ease of temperature control. 


Because the dielectric constant of biological materials 
even varies for samples from the same piece of tissue, 
accuracy and precision are not of primary importance. 


Furthermore, because the main concern in this work is with 
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microwave heating, and because much of this heating is done 
in the ISM (Industrial, Scientific and Medical) band at 2450 
Mhz, the frequency range from 2 to 3 Ghz is of special 


importance. 


What is important in biological measurements, however, 
is sample size and configuration, and ease of temperature 
controi. The method was therefore selected with these 


particular requirements in mind. 


In general, standing wave methods require carefully 
Shaped samples of rather large volume and accurately known 
dimensions, because the samples are required to fat 
appropriate waveguide structures. At the frequencies of 
interest, the apparatus would be bulky and difficult to 
incorporate into a suitable temperature control system. On 
the other hand, time domain techniques, that of Iskander and 
Stuchly (35) being a good example, require only small 
samples and could easily be incorporated into a temperature 
control system. The disadvantage of this technique is that a 
Sampling oscilloscope and an appropriate facility for 
Fourier Analysis are required. A time domain reflectometer 
with good spectral density at frequencies to 3 Ghz is also 
necessary. Because these facilities were not available, the 
cavity perturbation technique was chosen. This method 


adequately fulfils above criteria. 
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Of the many perturbation techniques available, three in 
particular were examined. The first of these was developed 
by Rueggenberg (36); the samples under study need not 
necessarily be smail as usually demanded by perturbation 
theory. Furthermore, samples of arbitrary shape can be 
measured, although calibration samples of identical shape 
and known dielectric properties are required. This method is 
potentially useful for measurements of whole organs. The 
second method is due to Rzepecka (37); its chief advantage 
is that the cavity employed is a section of waveguide, with 
attendant ease of fabrication. Special arrangements for 
temperature control in this method are also a desirable 
feature, but measurements on high-loss materials do reguire 
very small sampies. The third method examined, and the one 
finally chosen, was that of Risman and Bengtsson (38), 


described below. 
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The principal advantages of the cavity perturbation 
method developed by Risman are enumerated below: 

(1) The cavity employed is reasonably small and can 
easily be placed inside a temperature controlled 
environment, permitting control of sample 
temperature. Coaxial cable feed to the cavity helps 
in this respect. 


(2) The samples themselves are contained in glass or 
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teflon tubes which hold about 2ml of liquid. Coring 
techniques could easily be used for solid samples, 
as described by Bengtsson (39). The use of teflon 
tubes is necessary when freezing water, in order to 
prevent tube breakage. 

Careful cavity design minimises the existence of 
spurious resonances, which reduces the possibility 
of making incorrect measurements. 

The measurements themselves are based on a series of 
easily obtained calibration liquids, thus ensuring 
<+5% accuracy throughout the range of €* and €". 

In addition to the cavity, oniy a variable frequency 
source, some means of frequency measurement (good to 
about 1Mhz at 2450 Mhz) and some means of measuring 
cavity attenuation are required. Absolute frequency 
measurements are not necessary since it is the 
frequency shift which is important, and so, _ few 
requirements are placed on the frequency measuring 
equipment. In this particular case, use was made of 
the sweep-reference output of the source. This 
output provides a voltage which is proportional to 
frequency, and it proved possible to use it, rather 
than a more direct method of frequency measurement. 
Actual measurements are very easy to make, 
consisting of an adjustment for resonance followed 
by a frequency and attenuation measurement. 


Simple graphical techniques could be used to obtain 
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the values of €* and €" from the measurements of 
frequency shift and attenuation, though these would 
be tedious for the many measurements contemplated. 
However, computer technigues were used here. 

(8) A cavity designed and built in Sweden by Risman and 
his associates was available. This reduced local 


manufacture of a cavity to a simple copying process. 


(2.3) Theoretical Description of the Method. 

The cavity chosen was cylindrical of diameter 96 mm and 
height 178 mm, and was operated in the TM012 mode. Samples 
of up to 5mm in diameter were allowed. The diameter to 
height ratio was chosen to avoid other resonances. For the 
Same reason, vertical slots were made in the cavity walis 
and the cavity was divided into three sections, each section 
being electrically insulated from the others. Insulating 
discs anda teflon cone axially supported the sample within 
the cavity. The cone automatically centred the sample within 
the cavity. The inner surfaces of the cavity were gold- 
plated. A removable top facilitated sample insertion. The 
physical arrangement of the cavity is shown in figure 2.1. 
Figure 2.2 shows the field configuration within the cavity. 
The E-field is parallel to the sample and vanishes at its 
ends. Because the tangential component of the E-field is 


continuous at the boundary between the air and the sample, 
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there is good field penetration into the sample. 


The coupling loops are arranged as shown in figure 2.1. 
Investigations have indicated that coupling varies with €! 
in the most favourable way, since cavity transmission 
increases with €* for fixed €" of the sample. This gives 
good resolution for a wide range of €* because €" is often 
low when €* is low. Cavity perturbation methods are, in 
practice, not absolute methods, and so calibration 
procedures must generally be used. Using a number of weil 
defined reference substances, Risman (38) established the 
following relationships between €', €", the frequency shift 


and the attenuation. 
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where F, and F, are functions of the indicated variables. 


According to first order perturbation theory, F, Seed 


constant and the factor (at) °*° 


can be derived assuming 
weak coupling and Q(unloaded) >> Q(loaded). Typical curves 
obtained using these assumptions are shown in figures 2.3 


and 2.2. Details regarding the liquids’ used in the 


calibration procedure are covered in Chapter 3. 


The accuracy of the method is hard to estimate, but 
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measurementS on water over a wide temperature range 
correlates better than +5% with von Hippel's data (40) when 
corrected for frequency (see figures 4.1 and 4.2 in Chapter 
4). For other liquids, especially those with the higher 
values of €* and €", inaccuracy is about +10%, becoming 
about +20% at points far removed from calibration points. 
The precision or reproducibility of the method is good, the 
spread being on the order of 1-2% for liquids at the same 


temperature. 
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Chapter _3. Dielectric Measurement Technique 


(3.1) Description of Apparatus. 


Figure 3.1 is a block diagram of the measuring 


equipment. The cavity is the one discussed in Chapter 2. 


A signal at the resonant frequency of the cavity 
containing the sample (2500-3000 Mhz) is fed through the 
isolator and cavity and detected on a power meter. The use 
of a power meter is preferred over the use of a crystal 
detector because linearity is easier to achieve over the 
wide dynamic range required. A voltage proportional to 
frequency is fed from the Signal generator to a 
potentiometer circuit as shown in figure 3.1. When the 
voltage derived from the reference voltage and the 10-turn 
pot is equal to that derived from the Signal generator, the 
reading on the meter is reduced to zero. At any other point, 
the meter reads full scale. Tests using a frequency counter 
confirmed that this circuit is consistently capable of 
giving frequency shift readings to an accuracy of better 
than 1 Mhz, which is more than adequate for the dielectric 
measurements. Figure 3.2 is a plot of the reading on the 10- 
turn pot vs. the output frequency of the generator as 
measured on a frequency counter. Because the source is a 
voltage-tuned backward wave oscillator, any fluctuations in 


the sweep reference voltage will also appear as changes in 
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frequency. Therefore, long-term fluctuations in measurements 
of frequency shift by use of the sweep reference voltage are 


expected to be small, and this is indeed the case. 


The cavity and its coaxial air-line feed cables are 
contained in a Delta-Design temperature controlled oven. Air 
line is used here because of its ability to withstand wide 
Variations in temperature. All R.F. interconnections are 
kept as short as possible to reduce losses to a minimum. The 
operating temperature of the oven can be controlled to 
within +1°C over the temperature range -80° to +500°C. 
Liquid carbon dioxide is used to achieve those temperatures 
where cooling is necessary. Heating is by means of an 


electric element. 


(3.2) Temperature Measurement _and Control. 

Sample temperature is determined by a temperature- 
controlled oven. However, because the sample is contained 
within a brass cavity with rather restricted airflow, it is 
found that up to an hour is required for the sample 
temperature to reach the temperature of the oven, a time 
which was unacceptably long considering the number of 
measurements that were to be made. Also, the long 
temperature equilibration time would require large 
Gfanes mies of expensive coolant for temperatures below room 


temperature, and would result in excessive evaporation at 
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the higher temperatures. In addition, measurements in the 
vicinity of the freezing point of a mixture of liquids are 
difficult because, in spite of constant temperature, the 
dielectric constant changes slowly. This phenomenon is 
thought to be due to ice-crystal formation causing changes 


in the proportions of frozen and liquid materials. 


Because of the above problems, a dynamic method is 
required whereby constant cooling or heating are maintained 
while Making measurements of dielectric constant at 
appropriate temperatures. Unfortunately it is impossible to 
make simultaneous temperature and dielectric measurements 
because of the perturbing effect of the temperature 
measuring device. It is therefore necessary to use the 
following steps to obtain the required data: (a) A cooling 
and heating procedure which gives the desired range of 
temperatures is established. (b) Using this procedure, 
measurements of temperature vs. time are made. (c) Using an 
identical procedure, measurements of dielectric constant vs. 
time are made. (d) Removing time from the two sets of 
measurements, curves of dielectric constant vs. temperature 


are plotted. 


Although three replications of temperature vs. time 
measurements, and six replications of dielectric constant 
vs. time measurements are performed, for each sample, only 


about half the time is reguired for this technique as 
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compared to the point-by-point method. The repeatability of 
the temperature vs. time measurements indicates that the 
heating and cooling rates for given oven settings are within 


+3%, thus justifying this temperature measurement procedure. 


Separate temperature vs. time curves for each sample 
are necessary because the specific heat, latent heat and 
freezing point of different samples varies considerably. 
This is only so in the range -60° to +23°C. Tests revealed 
that for the higher temperatures, differences are less than 
+5% and only one set of temperature vs. time curves are used 


for all samples. 


Temperature measurements are performed by passing a 
thermocouple through one of the slots in the cavity and 
immersing it halfway into the sample. For measurements below 
room temperature, the oven is set at -80°C and temperature 
readings are taken every minute for 11 minutes as the sample 
cools down from room temperature. After 11 minutes, the oven 
is reset to #40°C and readings are again taken every minute 
for a further 23 minutes. These times are chosen to give a 
binimum sample temperature of about -60°C and a maximum of 
about +30°C. Equilibration at the low temperature is not 
attempted in order to conserve coolant. For measurements 
above room temperature, readings are taken every minute for 


16 minutes with an oven setting of +200°C. 
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Identical procedures are adopted during the dielectric 
measurements, with thermocouple removed from the sample. 
Cross-correlation of the temperature and dielectric data is 
achieved by the use of forms of the type shown in table 3.1. 
Except for the region around the freezing point, where 
different degrees of super-cooling sometimes occurred, 


overall temperature measurement accuracy is +5%. 
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In order to calibrate the cavity, a number of standard 


solutions were made (38). These solutions and their 
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dielectric properties are detailed in table 3.1. 


Using these values and the measured values of 
attenuation and frequency shift, as read on the power meter 
and the 10-turn potentiometer, a set of calculation graphs 
are constructed (38). Figure 3.3 shows €' vs. frequency 
shift, and is constructed by simply plotting the appropriate 
values. Figure 3.4 and 3.5 are constructed in conjunction 
with one another. Measured values of attenuation are first 
plotted on figure 3.5 along with their corresponding values 
of €". A reference line of slope m=(10g100-10g10) /(30- 
10)=0.05 (all logs are to the base 10) (see equation 2.2) is 
drawn for each plotted point in figure 3.5, and its 
equivalent value at an attenuation of 25 dB (i.e. €"(25)) is 
found. This value is transferred to figure 3.4 where it is 
plotted against log(€'-1). Using the values of €"(25) for 
any desired €", the other sloping lines of figure 3.5 can 
now be drawn parallel to the reference line. Graphical 
determination of €* and €" from measured values of frequency 


shift and attenuation then comprises the following steps: 


(1) Using the freguency shift, enter figure 3.3 
and obtain €'. 

(2) With this value of €", the corresponding €"(25) 
is found from figure 3.4. 

(3) Using €"(25) and the measured attenuation, €" 


is then found using figure 3.3. 
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This procedure is very tedious in view of the fact that 
these operations have to be repeated three times for each 
point to be plotted in order to give values for the mean and 
standard deviation. To facilitate the calculations for the 
many points required, computer programs to perform the above 
Operations were therefore written. These are described in 


detail in section 3.9. 


(3.4) Corrections for cavity expansion and contraction. 


The resonant frequency of a cavity is primarily 
dependent on its dimensions, and, being made of metal, the 
cavity dimensions are subject to variation with changes of 
temperature. In a situation where the parameter being 
measured is frequency shift, the temperature at which the 
loaded and unloaded resonant frequency is measured is not 
important, provided it is the same in both cases. However, 
in the present application, situations arise where the 
unloaded measurements are made at room temperature, while 
the loaded measurements are made at temperatures in the 
range -60°C to +100°C. To compensate for the resulting 
dimensional changes, appropriate corrections as shown in 
figure 3.6 are applied. This curve is obtained by measuring 
unloaded frequency as read on the ten-turn pot while 
subjecting the cavity to the heating and cooling procedure 
of section 3.2. Because sample and cavity temperature 


differ, cooling or heating time is again used as the common 


i Nae 


" 7 ae erus 7 


e ie eh ott Ro ae r " e198) 
yw <o% await gout Sets beat va od 9 
iN 


a 


, 


‘ma geen ods t02 saab ; vig 4 2 pa 


oO 
1 
Las 


add 20% anoitaius ino eat arBst red 


gyows G42 wa tISy OF ame 2eU34 slag ysiieae wn 
ai sdesocel Wes <Gepdal .fetsEae oo eists #50" oaenees 


7 ae songooa : 


| ri 
etrusetay @2 o¥aive>a « io yousopaat FORLORUE 
af? .tagon. to shee eaied ,bns lean cdi ask: 7 ) Hite aah 
30 A sgnants daky (apifaitsy of OST eeeS ae svete (3069, 
paded tofeaeian 949 8ISAs Hod douse. B m2 2 <0 
eit: Aoidw G6 S25 1679909 aa? “res az tonewpes’ | 4 at pea 
eon ai hengesen =f (paetpeTs +ngnogan conten aa re 
~aeveuoll .2e@a> A700 i ose sit 2k 5}: babavory ‘wha satin 


te 


eit asodw seins aaotzeutis aoksanthngs titeas som 
Z 

a 

elidy ,otsPatagav? §003 Te shea 548 atheessumuom ) a hee so r 


edt 46 aaqusasagiey (6 show? 83 meres bo bae sik re 
pate lvaces a7 FOR siscosqgos OF segure. 7% 0 

all uwote 2% ‘eerz2v91T0D «6|te 

; 

etsy #04 auton ous ao 

be eos ana pelsaed ada ine 


- nee 5 -. wn a - 


oe 


element to correlate correction values with measurements. 
Though sample temperatures at these times are different for 
different solutions, the effect of the samples on cavity 
temperature is negligible, and the same correction curve 


applies to all measurements. 


There are three different experimental procedures. 


These are: 


(1) Varying temperature measurements from room temperature 
to -60°C. 

(2) Varying temperature measurements from room temperature 
to +100°C. 

(3) Varying concentration measurements at several fixed 


temperatures. 


The experimental procedure for each type is slightly 


different, and will be described below. 


(3.6) Varying temperature measurements (-60°C to +20°C). 


For each different solution, three temperature 
calibration runs are performed in the manner described in 


section 3.2. Teflon sample tubes are used for distilled 
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water and 1% DMSO, and quartz glass tubes are used in ail 
other cases. Each run uses a different, newly mixed, sample. 
The mean temperatures obtained are entered in the 
appropriate column of a table like that shown in table 3.1. 
Appropriate correction values for frequency shift are also 
entered in the table. With the thermocouple removed from the 
cavity, and with an empty sample tube in place, the signal 
source is adjusted for maximum transmission through the 
cavity. The signal level is then adjusted to give a 
convenient reading on the power meter (usually +10dBm), and 
the zero set control on the frequency deviation meter is 
adjusted for balance with the 10-turn pot set to zero. With 
a newly mixed sample in place, the 11 minute cooling and 26 
Minute heating cycle is started. At one-minute intervals, 
the Signal source frequency is adjusted for maximun 
transmission and the values of frequency shift and 
attenuation are recorded on the table. To ensure readings 
exactly on the minute, frequency adjustment is started 
Shortly before the minute, with the reading being taken at 
the appropriate moment. Because rates of change are slow, 
the error attributable to this source is negligible. The 
procedure is repeated for 6 different newliy-mixed sampies. 
The samples are newly mixed to reduce the possibility of 
evaporation, which would cause an increase in the DMSO 
édnGbn tration: the latter being less volatile than water. 
Sets of 6 measurements are chosen aS a compromise between 


obtaining reasonable values of standard deviation and 
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available time. 


(3.7) Varying temperature measurements (+20°C_to_ +100°C). 


The measurements above room temperature are performed 
in the same way as in section 3.6, except that only one set 
of temperature vs. time measurements are made. Differences 
between the behaviour of different liquids is negligible as 
confirmed by temperature measurements on different samples. 
At these higher temperatures, some evaporation occurs, but 
in the relatively short measurement time involved, the 
expansion of the liguid, which in any case rules out the 
possibility of a cover, is more important. The overflow of 
the expanded liguid collected in the conical tube support, 
and it was thought that this might lead to inaccurate 
results. However, tests revealed that this is not the case, 
probably because the electric field at the ends of the tubes 
is very small. Similarly, removal of small amounts of liquid 
from the tubes did not affect the measured dielectric 
constant. The density change caused by liquid expansion 
always occurs for heated liquids so does not affect the 


results. 


Varying concentration measurements at room temperature 
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(23°C) only require obtaining balance with the ten-turn pot 
with and without a sample. In the case of the 5°C and 50°C 
measurements, however, a set of 11 tubes containing 
solutions of 0-100% in 10% steps are placed in the oven 
along with the cavity. After an equilibration time of about 
1 hour at the desired temperature, the equipment is balanced 
in the usual manner, the oven opened, and the empty tube 
substituted for one of the 11 samples. After a further 10 to 
15 minutes eguiiibration time, the appropriate measurements 
are made, and then the process is repeated for the second 
Sample. After one measurement on each sample, the equipment 
is again balanced, and a second measurement made on each 
sample. At this point, a new set of samples is prepared, and 
the process repeated. In this way, six measurements are 
obtained at each percentage concentration. In the case of 
the 50°C measurements, the tubes are sealed with tape to 
reduce evaporation during the four hour period in the oven. 
To allow for expansion, the tubes were not completely filled 
at room temperature. There was no observed evaporation at 
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(3-9) Data Reduction. 


Computer programs were written in FORTRAN for data 
processing; these are listed in appendix A and only briefly 
described here. Three programs are used in the data 


reduction process. Two programs are for data analysis and 
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one is a plotting program which graphs the values of €*', eE" 
and loss tangent vs. temperature. The plotting routine is a 


simple modification of a library program. 


The two analytical programs are a calibration routine 
and a routine to convert the measured values of frequency 
shift and attenuation to values of €* and €". Their function 


is described below. 


(1) Calibration routine: This program is called CALCOEFF, 
its purpose being to calculate the regression coefficients 
for the curves of figures 3.3 and 3.4. In order to achieve 
this, the program takes the measured values of frequency 
shift and attenuation, six of each in this case, and 
calculates a mean, an upperbound and a lowerbound. Combining 
these values with the known calibration values of €* and e€e" 
results in 3 sets of coefficients for the curves of €' vs. 
frequency shift and another 3 sets for the curves of log (€'- 
1) vs. Logé"(25). The coefficients are obtained using a 
bicubic spline fit of the measured data. Figures 3.7 and 3.8 
compare the measured data with the curves calculated from 
the spline fit for each of the two cases. These regression 
coefficients apply to the frequency shift data as obtained 
with the measuring system of figure 3.1. They can be used 
with any frequency measuring system by relating data from 
such a system to figure 3.2, although this would only apply 


to the particular cavity used here. 
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(2) Calculation of dielectric values: This program is called 
CALCULATE and its purpose is to obtain values of €' and €" 
from measurements of frequency shift and attenuation on 
unknown liquids. Aithough the temperature and concentration 
at which the measurements are made does not enter directly 
into the calculations, these values are also carried by the 


program to enable the appropriate printouts to be made. 


The program CALCULATE executes the following steps: 
(a) Read the temperature or concentration with the 
corresponding freguency shift and attenuation. In this case, 
6 values are read each time. 
(b) Using these values and the three sets of regression 
coefficients for figure 3.7, three values of €! are 
calculated at each value of temperature or concentration. 
These values represent the mean and upper and lower bounds 
of €* for that point. 
(c) Using the calculated values of €*, and the regression 
coefficients of figure 3.8, three values of E€" are obtained. 
Again, these values represent the mean, upper and lower 
bounds. 
(d) The mean, upper and lower bounds of €" are now 


calculated from: 


€" = alog(AT/20+loge" (25)-1.25) eA as.1) 
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where AT is the attenuation. This equation is a mathematical 
representation of figure 3.3. 
(e) For each temperature or concentration, print out the 


three values of €' and the three values of €*, 


The program is written so that it will accommodate any 
specified number of different temperatures Or 
concentrations, and any specified number (two or more) of 
repeated measurements. The accuracy of the algorithm was 
checked by obtaining values for selected points using the 
graphical technique, and comparing these values with the 
corres ponding computed values. There was never more than 1- 
2% variation, which is easily attributable to the 


inaccuracies of the graphical technigue and the spline fit. 
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(4.1) Presentation 


In the course of the dielectric measurements, a 
considerable amount of data was obtained and a brief 


discussion concerning data presentation follows. 


Curves of e€' and €" are presented for distilled water 
and for 1%, 5%, 10%, 20%, 70% and 100% DMSO. These comprise 
figures 4.1 to 4.14. All the curves in this group are drawn 
with the same axes to facilitate comparison. Parts of the 
curves for distilled water are shown dotted on those for 
each of the other solutions to further help comparisons of 


the behaviour of the different liquids. 


In each case, the error bars represent one standard 
deviation as computed on the basis of 6 separately measured 
values. Some error bars have been omitted, particularly 
those in the vicinity of the freezing point. The reason for 
this is that the rapid changes encountered in this region 
cause small uncertainties in temperature to become large 
uncertainties in dielectric properties, resulting in large 
and confusing error bars. In addition, super-cooling effects 
result in erratic and unpredictable behaviour in the 


freezing zone. 
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A further problem encountered was that the curves for 
dielectric properties obtained during cooling differed 
somewhat from those obtained during heating. The reason for 
this is that, as the freezing point is approached fron 
temperatures either above or below, reasonable homogeneity 
is maintained in the sample. However, after passing the 
freezing point, a mixture of solid and liquid is established 
because the liquid does not instantaneously freeze 
throughout, due to latent heat requirements. The dielectric 
properties of this mixture of phases are ill-defined and 
certainly different to those of the homogenous material. 
Therefore, to avoid confusion on the graphs, only points 
plotted (i.e. points above the freezing point during 
cooling, and points below the freezing point during 
heating). In some instances where super-cooling is involved, 
this policy necessitated interpolation between adjacent 
values in order to obtain a value right at the freezing 
point for the cooling curve. This point is then connected to 
the value for the same temperature as obtained during 
heating. Errors_in the immediate vicinity of the phase 
change are of the order of 450%, and caution must _be 
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The second group of curves shows the behaviour of the 
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dielectric properties of water and a 5% (w/w) solution of 
dextran. These comprise figures 4.15 and 4.14. The axes in 
this case are different to figures 4.1 to 4.14, hence the 
repeated curves for water as a standard of comparison. A 
maximum temperature of 40°C is used here because dextran is 
primarily used in applications at temperatures below this 
value. The comments made above with regard to the first 


group also apply. 


The third group of curves comprise dielectric values 
for 10% DMSO and solutions of 10% DMSO with 0.1% Nacl, 1.0% 
NaCl and 5% dextran (w/w). These are figures 4.17 and 4.18. 
Again, the axes are different from those of the first group, 
and so 10% DMSO is repeated for comparison purposes. 
Previous comments with regard to maximum temperature also 


apply here. 


The fourth and final group of curves, figures 4.19- 
ee comprises measurements hade for different 
concentrations of DMSO at various fixed temperatures. In 
this instance, values of the loss-tangent are also plotted 
for reasons made clear in section 4.2. Loss-tangent values 


were not considered necessary in the first three groups. 


(4.2) Comparison with Literature for Distilled Water. 


In order to ascertain the accuracy of the measurement 
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technique, and also for purposes of comparison, dielectric 
measurements on distilled water in the temperature range 
from -60°C to +100°C were made. The results of these 
measurements along with those of von Hippel (40) are 
presented in figures 4.1 and 4.2. For both €* and €", +1% 
agreement is indicated except at high temperatures for €! 
and in the immediate vicinity of the freezing point for €", 
where the deviation is approximately +10%. This suggests 
that the technique employed in the present work provides 
results which are consistent with those in the literature. 
In this context, it should be pointed out that for very low 
loss materials such as ice from distilled water, the values 
of €" are not expected to be accurate unless the calibration 
curves are adjusted for the presence of the sample holding 
tube. The values presented in appendix D for ice are 
therefore not accurate, but represent the minimum resolution 
of the calibration curves as obtained in Chapter 3. The 
dielectric loss for ice is too small to appear in figure 
4.2, and in the case of the other mixtures, losses in this 
region are higher, and an accuracy of better than +10% is 
expected. In addition, very high loss materials (loss 
tangent>3) introduce considerable attenuation and this leads 
to measurement difficulties because the attenuation readings 
are close ys the noise level, and accurate tuning for cavity 
resonance is difficult. This explains the larger error in 


the vicinity of the freezing point. 
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Moving now to the results obtained for the first group 
of measurements, it is seen that the behaviour above the 
freezing point of €' for all solutions except those of 70% 
and 100% DMSO is very similar to that of water. The freezing 
point for any DMSO solution can be established by reference 
to the phase diagram of figure 4.22 (41). In the case of €", 
the behaviour is again similar but with generally higher 
losses for the DMSO solutions. For all solutions excepting 
70% DMSO, the values of both €* and €" decrease with 
increasing temperature above the freezing point. In the case 
of €' for 70% and 100% DMSO, there is a considerable 
Variation from the corresponding values for water, although 
the behaviour with temperature is similar. €" for 70% DMSO 
shows a broad peak around 23°C, a fact which will be 
discussed in more detail below. For the solutions on which 
measurements were made, the maximum Loss occurred in the 
case of 20% DMSO at a temperature just above the freezing 
point. This can be explained by consideration of the fact 
that liquids have what might be considered an instantaneous 
structure (42) comprising, in the present case, hydrogen 
bonded molecules. The proportion of bonded molecules at any 
instant increases as the temperature is lowered towards the 
freezing point, and so the molecular forces impeding the 
dipole orientation dominate, and the dipoles become less 


able to follow the changes in applied electric field. Hence, 
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an increasing phase-lag between the field and the dipole 


orientation develops, and losses (i.e. €") increase. 


In all cases, except for 70% DMSO, a fairly abrupt 
change in both €* and €" is evident at the freezing point. 
This may also be the case for 70% DMSO, but the freezing 
temperature of -72°C was below the measurement range. These 
abrupt changes are associated with the phase change at which 
the continuously changing instantaneous structure of the 
liguid becomes a rigid solid. The effect of this is to 
reduce the relaxation frequency. For instance, the Cole-Cole 
plot for ice (43) (figure 4.23) shows that the relaxation 
frequency of this material is about 3.0 Khz as compared to 
18 Ghz for water (figure 4.24) (43). At frequencies around 3 
Ghz, as used in the present measurements, €" for ice is very 
small and €'3.2 and _ shows very little temperature 
dependence for the range of temperatures considered here. In 
the case of the DMSO mixtures, however, losses are 
considerably increased over those of ice, even in the case 
of 1% DMSO. In addition, the change in the dielectric 
properties as the freezing point is approached from below is 
relatively gradual. These phenomena suggest that hydrogen 
bonding between the water and the DMSO is taking place, 
thereby reducing the relaxation frequency of the bound 
water, and therefore increasing €". This also occurs in 
other situations, for instance, where water molecules form 


bonds to the cellulose in wood (44). The cryobiological 
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implications of these results will be discussed in section 
4.3. 


(4.4) Comparison of Cole-Cole Plots for Water_and_ DMSO. 


The fact that €" decreases with increasing temperature 
above the freezing point suggests that the measuring 
frequency is below the relaxation frequency for all the 
concentrations considered in this group, except for 70% DMSO 
(figure 4.12). For the liguid phase, this can be explained 
by reference to figure 4.24 which is a Cole-Cole plot for 
water and for DMSO drawn on the basis of literature data 
(16,17,43), and some of the present measurements. It is 
evident from this figure that the relaxation frequency for 
water is in the vicinity of 18 Ghz, and that for DMSO is in 
the vicinity of 9 Ghz. For measurements at frequencies below 
the relaxation frequency, decreases in both €* and €" with 
increasing temperature would therefore be expected. This 


corresponds to moving from curve A to curve B. 


This point is further clarified by reference to figure 
4.25 (40). Here it can be seen that the relaxation frequency 
increases with increasing temperature. Considering a 
measurement at frequency f, below the relaxation frequency, 
it is clear that a drop in value from €",, to €",2 would be 
expected for a temperature increase from T, to 1,2. On the 


other hand, for a measurement well above the relaxation 
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frequency, such as fag, an increase from €"gq, to EM", is 
expected for an increase in temperature. From this 
discussion, it is apparent that a peak in the value of €" as 
the temperature is varied would suggest that the measuring 
frequency was in fact the relaxation frequency at the 


temperature of the peak. 


(4.5) Relaxation in _a_70% Solution of DMSO in Water. 

Further consideration of figure 4.24 suggests that 
increasing the DMSO concentration up to 70% is analagous to 
increasing the measurement frequency or reducing the 
relaxation frequency of the solution. Increasing the 
concentration from 70% to 100% has the same effect as 
increasing temperature, and from the above discussion, this 


has the effect of increasing the relaxation frequency of the 
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This implies that the hydrogen bonding between the DMSO 


and the water molecules reaches _a_maximum for the eutectic 
mixture (i.e. approximately 70% DMSO). Szmant (14) has noted 
several other effects that occur for the eutectic mixture, 
and he indicates that the addition of DMSO to water 


increases the organization of water up to concentrations of 


0.33 mole fraction DMSO (i.e. 66%). The composition of the 
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liquid is then DMSO.2H,0. 


Referring to figure 4.12, a definite peak in the 
temperature curve of €" occurs in the vicinity of 23°C. This 
leads to the conclusion that the relaxation frequency for 
70% DMSO at 23°C is at the frequency at which this 
measurement was made (i.e. 2750 Mhz). Because of the large 
error-bars in the vicinity of the peak (due to a large 
attenuation through the cavity), the exact position of the 
peak is ili-defined. In addition, the freezing point of 100% 
DMSO is at 18°C, only 5°C below this temperature, and so the 
values at temperatures below the maximum may be due to some 
partial freezing phenomenon in the liguid. The even larger 
error-bars in this region tend to support this possibility 
because partially frozen solutions are characteristically 
ill-defined as to their composition. On the other hand, the 
fact that €" for other concentrations continues to increase 
right down to the freezing point reinforces the conclusion 
that the 70% solution does exhibit a minimum relaxation 


frequency. 


(4.6) 5&2 Dextran Solution. 


Figures 4.15 and 4.16 are plots of the temperature 
behaviour of a 5% (w/w) solution of dextran. Dextran is one 
of a class of high molecular weight cryoprotective agents. 


The results show that the effects of dextran on the 
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dielectric constant are minimal, although the solution is 
slightly more lossy than water above the freezing point, and 
the transition beiow the freezing point is not quite as 
abrupt. Below this transition region, the solution behaves 
very similarly to pure water, suggesting that the 
macromolecules have little effect on the ice structure. 
Because of the large size of the molecules, it is likely 
that more significant effects would be apparent at lower 
frequencies, as has been indicated on proteins and other 


biomolecules (43). 


(4.7) 10% DMSO_and Various Additives. 

Figures 4.17 and 4.18 are the results of measurements 
on solutions of 10% DMSO and mixtures of 10% DMSO with 
either 0.1% NaCl, 1.0% NaCl or 5% dextran (w/w). Most 
biological solutions contain ionic components and so the 
NaCl was added to give some indication of the effect of an 
ionic component on the dielectric properties. As can be seen 
from the results, the only significant effect is the 
addition of 1% NaCl to the 10% DMSO, and it is only €" above 
the freezing point that is affected in this case. 1% NaCl is 
close to the widely used physiological saline solution of 
0.9% NaCl, and was chosen for this reason. The increase in 
conductivity resulting from the addition of the NaCl 


accounts for the increase in €" as observed in figure 4.18. 
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(4.8) Behaviour of €" with Concentration. 


The final group of curves comprises figures 4.19, 4.20 
and 4.21. In these, the dielectric properties, including the 
loss tangent, are plotted against the DMSO concentration for 
three different temperatures. Figure 4.19 shows that the 
behavior of €" at 23°C and 50°C is as one would expect on 
the basis of linear mixing, that is, with the values 
decreasing linearly from those of water to those of pure 
DMSO. The curve for 5°C (figure 4.19), however, has a 
minimum value at 70% DMSO. This is consistent with the curve 
of figure 4.11. The abrupt change at about 90% DMSO is due 
to higher concentrations having freezing points above 59°C. 
The reason for the minimum in the value of €' at 5°C can be 
explained by reference to figure 4.26, which is a typical 
Debye curve of the behaviour of a dielectric. Previous 
results indicated that DMSO at a concentration of 70% has a 
minimum relaxation frequency, fp. Therefore, at 10%, 
measurements are being made in the vicinity of point f7 on 
figure 4.26, with a resultant value of €"5. For any other 
concentration, however, the relaxation frequency is higher 
which is analagous to being in the region to the left of f, 
such that f<f>. For any point in this region, the value of 
€' will be greater than €'> and so one would expect €'> to 
be the minimum value. The fact that this is so for all the 


concentrations in figure 4.19, except those below the 
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Turning now to figure 4.20, we see that €" goes through 
a maximum at about 40% DMSO, for a temperature of 5°C, but 
has maxima at 70% for 23°C and at about 50% for 50°C. 
Reference to appendix D shows that the high losses in the 
vicinity of the maxima of these curves results in large 
Standard deviations and so the exact locations of the maxima 
are hard to determine. However, where comparison is possible 
with the earlier curves, the data is consistent. The trend 
towards increasing losses with decreasing temperature and 
increasing DMSO concentration iS maintained throughout 
except in the vicinity of 70% DMSO where the behaviour is 
consistent with a dispersion region. Figure 4.21 is a plot 
of the loss tangent vs_ concentration, and for all 
temperatures plotted, it shows maxima in the vicinity of 70% 
DMSO. This can be explained by considering the definition of 
loss tangent (e"7e'). Figure 4.26 shows that E€"/e" could go 
through a maximum in the vicinity of f,, though this would 


depend on the behaviour of €* in this region. 


The results are now examined with a view to the effects 
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of microwave heating of biological materials containing the 
solutions as cryoprotective agents. As the results have 
indicated, the dielectric loss of frozen materials is 
considerably less than it would be at room temperature. This 
is generally true of all frozen materials containing water, 
and, specifically, of biological tissue. In cryobiological 
applications, the biomaterials to be heated are initially in 
the frozen state, and so the expected dielectric loss is 
small. This is in direct conflict with the requirements of 
rapid heating, and any means of increasing the dielectric 
loss at low temperatures is useful. From figures 4.3 to 
4.12, it is clear that the addition of DMSO increases the 
losses at low temperatures. In general, the higher the 


percentage of DMSO, the higher the increase in the losses. 


One of the important areas of current cryobiological 
research is in establishing the most effective 
cryoprotective agent. With regard to DMSO, various 
concentrations have been employed, usually in the range 1- 
10%. However, Pegg (45) has used apparatus in which the 
concentration is increased up to 70% as cooling proceeds. 
Aside from the cryoprotective properties, the question of 
toxicity arises, and concentrations above about 15-20% are 
usually considered to be toxic (15). However, reducing the 
temperature allows the use of higher concentrations, hence 


the technigue of Pegg (45) is acceptable. 


i 


ve a ie ale 


odd pakatssaod 2Lsbise en 
gysd 2tivasr ast 2h sada 
2h ‘elsizetna aexvowk 40 

sid? .a:eseuegesromoor 36 scl 


= 
Ki 
ek. 


ivter Cuintstgod sisite7am noe 


Lasiytiobioyio ot seezecd SH. pes i 2 


mi Yffskvink ore Dagped ot oF ata ver , 
ai enol Sbitselelhi fasvegus alt" ys ban wtee8 soxoat 99 ea ; 
to atnanosthipe2 odt atie Joattaes nik et afdt .ilsse -~ 


niygoaleib ad? pateeetons 10 a enna 
of Eo eotupe? sort .fetven contegacanes vol 36) aeol 
gue aeupetoni GEnd de sogstbis wkd s6ike zaohs si ti ashe 
sh} safghd of ,léceup at -eendatogied wol 46 contol: if 
~2ezanl ahs ai se nezodt omy yeepld’ aay ened iq epszas01eg | 1 


‘nes 
leakouloidayy> taeigvua 30 B5a56 snamzoges aids fo 220 | | 
avisoetie s200 ‘9h? antdetfasiee ai al do seal - 
auor psy orig o? Huse atli .7seps svitoes07gOYI9 y 
_} sugies add ab Ykeveo .b4yolque deed evéd eaodtea2ae0009" | 
a4 doidw di eetsaeqqe been ~ebd (28) eps? ,rsvewol kor 
.#bes002q gHileos es 0% of qu beasetomh et doi saxd.ag0G0> 
to sulstzadp dt ee itreqgoty svitastemger¢3o ods mozart obiek ie 
eis 200-2et tuadn evods ehottsrtagod6s bas ,eseias yoipixod 
add hiquher steeswow . (21) sixot Bd OF bereblanos {ihsoee 
enagd yatoitéatnene> iodeid To sal eat ewolis cient 

SP ldsiqaooe ei i tadal io sop hades sat 


; se » 2 ers aS 


7 a 


- 


ae | 4 
ne ea 


80 


Referring again to figures 4.3 to 4.12, it appears that 
maximum enhancement of iow temperature heating would occur 
with 70% DMSO. This is only practical when the concentration 
is increased with decreasing temperature, and so for the 
more usual application where a single concentration is used, 
a limiting concentration of 10% is indicated. Reasonable 
enhancement of the dielectric loss does occur at this 


concentration. 


As has been pointed out earlier, the critical range in 
tissue freezing and recovery is between -60°C and +10°C. 
Therefore, it is in this range that the possibility of rapid 
heating is most important, and so the low losses at 
temperatures below -60°C do not constitute a serious 
disadvantage. It may therefore be concluded that___the 
temperatures, and that higher concentrations of DMSO provide 


Another important question in microwave heating in 
cryobiology is that of heating uniformity, and the 
dielectric measurement results are now examined in this 
light. Consider a homogenous body which is at the same 
temperature throughout. The body is placed in a microwave 
field and heating is begun. It is inevitable that some non- 
uniformity of field distribution will exist, and this will 


cause different areas to heat at different rates. From the 
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results presented in the curves, and from other dielectric 
data on frozen biological materials (46), it is clear that 
the warmer the material, the higher the loss and therefore 
the more energy it will absorb. This will tend to increase 
the non-uniformity of temperature distribution, perhaps to 
the extent that temperature gradients are able to cause 
damage. This effect becomes more pronounced for those 
solutions in which the change in dielectric constant is more 
abrupt. Heat conduction effects are of course a moderating 


influence but non-uniformities will still arise. 


The dielectric properties of perfused biological 
materials are heavily dependent on the presence of the 
perfusate if this is DMSO (22). This means that the above 
behaviour will certainly be influenced by the properties of 
the DMSO solutions. It is reasonable to suppose that those 
solutions exhibiting the most abrupt changes in dielectric 
properties would cause the greatest non-uniformity of 
heating. In this category are water and the 1% solution of 
DMSO. Although not directly related to the present 
discusssion, the solution containing dextran also exhibits 
an abrupt change which suggests that cryoprotective agents 
of this type are less suitable than DMSO-based agents from 
the standpoint of microwave heating. On the other __hand, 


those__solutions containing higher percentages of DMSO_ (25%) 


show increasingly gradual transitions in the vicinity of the 


freezing point, and their use is__therefore suggested _from 
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In conclusion, then, DMSO solutions up to a 
concentration of about 70% (eutectic mixture), exhibit 
relaxation frequencies which fall from approximately 18 Ghz 
for pure water to about 3 Ghz for 70% DMSO. As _ the 
concentration is increased above 70%, the relaxation 
frequency increases to about 9 Ghz for pure DMSO. Of the 
solutions measured, the most lossy was a 40% solution at 
5°C, €" in this case being approximately 34. Of the other 
additives tested, only 1% NaCl had any appreciable effect, 
and then only on €" at temperatures above the freezing 


point. 


From a biolgical viewpoint, considerations of microwave 
heating rate and uniformity dictate the use of the highest 
possible concentration of DMSO permitted by considerations 
of cryoprotective action and toxicity. The technique of Pegg 
(45) where the DMSO concentration at low temperatures is 
high is certainly attractive from a heating standpoint. The 
implications of the results of the dielectric measurements 
on DMSO solutions underline the importance of a knowledge of 
the dielectric properties of cryoprotective agents when 


microwave heating is used. 
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The recent proliferation of high power microwave 
devices has aggravated the problems. of electromagnetic 
"pollution', but has also made new heating techniques 
available to the biomedical engineer. The effects of 
microwaves on living tissue have thus become a question of 
some importance, and the accurate determination of the 
electromagnetic fields within arbitrary-shaped inhomogeneous 
dielectric scatterers is a problem which is basic to this 
question. Broadly speaking, two approaches are available for 
the solution of this problem. These are: (a) direct or 
indirect measurement of the fields and (b) mathematical 
determination of the fields based on solutions to Maxwell's 


eguations. 


Some of the methods that have been used in the first 
category include: 
(1) The use of probes or loops to measure the 
electromagnetic fields at particular points. Because the 
presence of a metallic probe would considerably perturb the 
fields, non-metallic probes using liquid crystals to measure 
the heating effect of the field have been developed (47). 
Information is transferred from the probe to the measuring 


equipment using optic fibres to minimise field perturbation. 
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(2) The use of thin sheets of lossy material with which 
liquid crystals are in thermal contact. Colour photographs 
of the sheets after a period of heating can then be 
interpreted as electromagnetic field distributions (48). 

(3) Observation of the re-radiation from a diode implanted 
in the scatterer, which is in an electromagnetic field (49). 
(4) The use of colour thermovision on sectioned phanton 
scatterers to determine the heating distribution, from which 


the field distribution can be inferred (50,51). 


Though some of these methods have proved very 
effective, they are not as flexible as mathematical methods 
because of hardware requirements, and because of 
difficulties involved when measurements are to be made on 
biological tissues. For these reasons, it was decided to 
adopt a mathematical approach towards modelling the heating 


of frozen organs in a microwave field. 


In the past several years, a considerable number of 
mathematical field modelling techniques have appeared, and 
Bates (52) recently reviewed existing methods. For simple 
scatterers such as spheres, cylinders and cubes, analytical 
solutions which can often be obtained by separation of 
variables, are available. However, for the more complicated 
scatterers encountered in cryobiology, one must resort to a 
numerical method if an accurate model is to be constructed. 


Three classes of model exist. These are: 


Bis, Sone ae 
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(1) Transmission models: In these, assumptions about the 
reflection and refraction of the impinging waves are made, 
and the amount of energy absorbed in finite elements of the 
scatterer is calculated as the electromagnetic wave 
propagates towards the interior of the scatterer. The shape 
of the scatterer is limited to fairly simple cases, but 
different dielectric properties can be assigned to different 
regions. Typical of this type of model is that of Nykvist 
(28), which models the heating of a meat roast ina 
hicrowave oven, where energy is approaching from many 
different directions. 

(2) Frequency Domain Modeis: Most of the models appearing in 
the literature fall into this class, and the methods are 
based on the assumption of an exp(j2uft) time dependence in 
the fundamental Maxwell's equations. A set of linear 
eguations for either field variables or field expansion 
coefficients is derived, and then solved with a suitable 
matrix inversion scheme. Representative models of this type 
are those due to Shapiro (53) and Durney (54). In these 
models, the scatterers are restricted in shape and apply to 
a multi-layered sphere and a prolate spheroid respectively. 
The solution of a set of integral equations, as presented by 
Livesay (55) also falls into this category, but it is 
applicable to arbitrary-shaped bodies. The disadvantage of 
this method is that it would require the inversion of a very 
large matrix if it were applied to a large scatterer with 


fine structural details, as is the case in many biological 
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applications. 

(3) Time-domain Methods: These methods treat the irradiation 
of the scatterer as an initial value problem, considering 
propagation of waves from a source which is turned on at 
time t=0. A solution of a finite-difference analog of the 
time-dependent Maxwell's equations is obtained on a lattice 
of points including the scatterer. A good example of this 
technique was presented by Yee (56), being later extended to 
biological materials by Taflove (31,57). This method 
appeared to be the most suitable and was the one finally 
chosen for the heating studies performed in the present 


work. 


(5.2) Reasons_for Choice of the Method Used. 

For reasons outlined in Chapter 6, it was decided to 
investigate Microwave heating of canine kidneys under 
conditions of plane wave illumination. These organs have a 
complex structure with many fine details, and so the method 
chosen would have to be able to model these details with 
reasonable accuracy. In the model finally adopted, details 
down to Awm/20 were incorporated ()\m=Mininum wavelength). 
The possibility of approximating the kidney to various 
regular shapes such as spheres, cylinders and oblate 
spheroids was considered, but the resulting internal 
structure could not have had this degree of detail and still 


maintain the regularity and symmetry demanded by the models 
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intended for these shapes. The inclusion of the smail 
details in the model thus ruled out the first category and, 
excepting Livesay's model (55), the second category 
discussed above. Several factors influenced the choice of 
Taflove's model over that of Livesay. They are: 

(1) The accuracy of Livesay's model depends on the size of 
the scatterer in relation to the wavelength, which is not 
the case with Taflove's model. 

(2) To implement Livesay's mnodel, the scatterer would have 
to be sub-divided into N sub-volumes and three coefficients 
would have to be caiculated for each sub-volume. A numerical 
integration technique would probably be necessary for each 
of these calculations. The resulting 3Nx3N matrix would then 
have to be inverted to obtain the required solution. On the 
basis of data presented in Chapter 6, and assuming a maximum 
sub-volume size of one quarter of the minimum wavelength 
encountered, approximately 100 sub-volumes are required. 
Thus, a 300x300 matrix would have to be inverted. Although 
not prohibitive, this is a problem of some magnitude and 
would have no advantage over fTaflove's method in this 
regard. 

(3) A possible future extension of Taflove's method to non- 
plane wave irradiation of a scatterer surrounded by 
conducting walls (i.e. a microwave oven) appears to be a 
reasonable possibility (see Chapter 7). This would allow 


even more realistic modelling of microwave-heated materials. 
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(5-3) Description of the Method Chosen. 

As mentioned above, the method chosen is a time-domain 
approach which treats the irradiation of the scatterer as an 
initial value problem. At time t=0, a plane wave of 
frequency f is assumed to be turned on. The propagation of 
waves from this source is simulated by solving a finite 
difference analog of the time dependent Maxwell's equations 
on a lattice of points including the scatterer. Time 
stepping is continued until the sinusoidal steady-state is 
achieved at each point. The field envelope, or maxinun 
absolute value during the final half-wave cycle of time- 
stepping is taken as the magnitude of the phasor of the 
steady state field. The following description is based on 


that presented by Yee (56) and extended by Taflove (31,57). 


(5.4) The Yee Algoritha. 
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Maxwell*s equations in the rectangular co-ordinate system 


(X,Y,Z) can be represented as follows: 
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QEx 1 / Ys OH 
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DEG 1-H ea 
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> Fae) Hes Gah 


At this point, a set of finite difference equations 
following Yee's (Y¥3) notation are introduced. A _ space 
lattice point is denoted as: 

(i,jek) = (18, 355,kS) Pd hi) | 
and any function of space and time as: 
Foes kh) =k (Lord ag keene t) Ue beta oer (503) 
where =$x=sy=$z is the space increment and §$t is the tine 
increment. Yee uses finite difference expressions for the 
space and time derivatives that are both simply programmed 


and second order accurate in §& and $t. They are: 
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In order to achieve the accuracy of 5.4, and to realise 
all the space derivatives of 5.1a - 5.1f, the components of 
E and H are positioned about a unit cell of the lattice as 


shown in figure 5.1. To achieve the accuracy of 5.5, FE and i 
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are evaluated at alternate half-time steps. These 
assumptions lead to the following set of finite difference 


equations for 5.1a - 5.1f (56): 
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The assumption that ot/M(iejek)S is constant for all 


i,j,k of the lattice is now made, and the following 
constants are defined: 

R= St/E, ---.. S70 

Raz St? / (8°, é,) ae 6S Tb 

Rp = St / 4, & tress STE 

Calm) = to - Ro (m)/E,(m) «.---. 37d 


Cytm)= Ka /6-On) +... S726 


where m is a tissue type integer (simply an integer assigned 
to a given type of tissue). The proportional electric field 


vector is defined as: 


and so equations 5.7a -5.7e and 5.8 allow 5.6a to be re- 


written as follows: 
Here Ce jab Red) = Ho? (i, j+4, R42) + EY Cc, J+4, Ret) 


— Ey (je ) + Ey (4 Ret) — Ey (6, +l, R+ 2) sree ee Sg 


5-6b and 5.6c can be written in a similar way. In the same 


way, 5.6d can be written as follows: 
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5.6e and 5.6f can be written in a similar way. 


(5.5) Lattice Truncation Conditions. 


A basic problem with any finite difference solution of 
Maxwell*s equations is the treatment of the field vector 
components at the lattice truncation. Because of limited 
computer storage, the lattice cannot cover a sufficiently 
large portion of space so that the scattered wave at the 
lattice truncation might be closely approximated. 
Necessarily, the lattice must terminate close to the 
scatterer in a region where the nature of the scattered wave 


is unknown. 


Proper truncation of the lattice requires that any 
outgoing wave disappear at the lattice boundary without 
reflection during the continuous time-stepping of the 
algorithm. Improper truncation results in error for ali time 
steps after the boundary wave reflections return to the 
vicinity of the scatterer. One possible truncation condition 
is to set the field components at each truncation point to 
zero for all time steps. This situation is similar to that 
resulting if a conductor were present at that point, and 


would result in reflections. This is known as a “hard" 
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truncation condition. 


In the free-space situation being modelled here, the 
hard truncation condition is inadequate, and a "soft" 
truncation has been developed. One possibility for a soft 
truncation is to assign linearly increasing values of 
conductivity near the lattice boundaries, thereby absorbing 
the outgoing wave. For small reflection, however, this would 
require absorption over approximately 1 wavelength, which is 
intolerable from the point of view of computer storage. In 
addition, wavefront distortion for waves propagating 


parailei to the lattice boundary would be undesirable. 


To overcome this problem, a soft lattice truncation 
which relates the field values at the lattice boundaries to 
those 18 or 2§ within the boundary was developed. An exact 
condition is available for a one-dimensional lattice, since 
all waves would arrive perpendicularly at the lattice 
boundary. However, in the 2- or 3-dimensional case, one or 
several waves at differing angles of incidence could be 
arriving simultaneously, and no exact condition is available 
for this situation. Taflove (57) discusses this at length 
for the 1- and 2-dimensional cases, and only the 3- 


dimensional case will be dealt with here. 


A typical lattice is shown in figure 5.2. This lattice 


Was used to test the program on a sphere, as described in 


‘le tet ba) il. hae “- shaos:: sitate 
| eon | aa 
ait 10d tev t Hare pn tl at 6 vo ree pee 
sides, tance cama 
tion 5 ol yrShkdideog 400, .beq | aad gobs) | 
to contsy pataseapal gienesil & me md et andabomuad 
pelareads yiszads ,eeLsebened: a), tesa ysivisosbaos — 7 
Siwow alte daewed) wolt=eitve Linea A vovew gntonsue et 
et dotde .itpnelavey tT Livsee cco sya Gwe, Gontqnoads eakuyer ay 
vi .apnt0se retudeod LO ety to 14 | e032 oldsrelotat | 
paksopsqaty aevew aol er Samoa okt tb bs a 
shir besonw of SLido8-\Yaalhneod @odttut ad2 of foifsisq — +4 
ro 1 
404i )s0nUuTS SDEtABL Pins «4 vasidog arat smOsssvG. OT ©” i 
at ackrdbngod esitish oid ds devley Bett ods acseiss wnkse | 
snes 8a .beqolevel Sew ys6hspod ed% madgtw BS. 20> avian 
eonia ,eskdtel Ishoteechi-pw & 10% ekdgiiovws 24 noisibaoo — 
anittei att #6, yiaaly> isasg 16g eviaee biuow seven: fis” 
7060 co. geeso [endieasmih-f 340-4 odd at ie vewoll strebauod 
ot biueoo #oaebinal Yo ezelpia. , pat edhee +s «6eovew Leisven 
wtishiows me aod +2 thaws fasxe of bas’ santeans eee earn, TT 
itjynet *@ 9d Sepauvait (Ve) svelaa® sno ts0sbe, akid S05” 
~t off ghao Boe ,apecn Isnotenemlbepitay =f stent 
asd dtu cee ie snails ieaokanses6 | E 
; te Gay 


eolttal sid®-.¢.2 evveit al ovis ieevor ee me “i 


- 


section 5.9. From the basic time step relation of the 


algorithms: 
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That is, it takes two time steps for a wave in air to cross 


between two adjacent lattice points. The integer constant, 
Clim) = ZACH FZ ACM) eee GAD 


is defined as the number of time steps required to propagate 
between adjacent lattice points in tissue type ao. (1) is 
the wavelength in air, and (m) is the wavelength in tissue 


type m. Further, the stored field vectors 
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are defined. For the lattice plane x=1/2 (see figure 5.2), 


the truncation is given by: 
A A Sra 
Ho (a,j, Rrz)=DHy (2, ),R-2) + Hy Cz, J, Ree) 
+ Hy (3,5 R+2)4 Nise Ce ae saith IGS 


which in effect takes the mean of the H-field values at the 
point, and its two adjacent lattice points as existed 1 time 
steps before, and assigns that value to the lattice 
truncation point. Similarly, the other lattice truncation 
conditions are: 

For plane x=1/2 . 
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(5.6) The Symmetry Conditions. 


If even symmetry of the scatterer about one or two 
planes can be assumed, an important saving in computer time 
and memory space can be realised. In the case of the canine 
kidney, symmetry about two planes is assumed (see Chapter 
6). An incident plane wave having components E and 4H Ls 
assumed and so the following symmetry conditions can be 


defined. 
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The simulation of a continuous sinusoidal plane wave 
source is now considered. One possibility would be to vary 
the electric field along the points of one endface in the 
appropriate manner, but this would represent a "hard" 
truncation condition for any scattered waves, and so cause 
undesirable reflections. Alternatively, a sinusoidal 
Variation can be superimposed on the existing field at the 
points in a plane imbedded 2 or 3 into the lattice from 
one endface. This allows reflections to pass through the 
source and be truncated in the manner described above. This 
condition simulates a plane wave originating at infinity and 
a reflected wave returning to infinity without further 
interaction. The appropriate equation defining this type of 
source is: 

© (6,3, Reb) < loco Rp sin(2m4nSt) + EZ Le,8,Rrs) Ste 
This defines an impinging field with an E component being 


developed at plane 35. 


(5.8) Stability Criterion. 


Taflove (57) derives the stability criterion for the 


algorithm, and arrives at the following result: 


ai= 
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where Vmax +S the maximum phase velocity expected within the 


voN rea re Pea ae . * 


syew custq laiinepale avout 
e2sv of of OS haw retaadiesog ou ) 
oie ri wonthad én 30 atibeg ent 
“fHan/° ss isseeaged Bivoe 21nd? ns i a 
sapas ue Sab yeavew Bodgedeoe yrs 208 java Lhaeo pres 
iiicevers , & ylewkeentedts anpgsoatzes “pe 
ga 2% 5lske pnttaise S27 126 vacua tbanie ec S89 ‘aos ba us 
woot evicted bad osPE E se S bebhadait impel Bak _ atakog: ° se a 
o49 #paoddy Bese oF eoodynelies woe nidt oe py) 
aid? -9vods bednxsaw s9nTse ed? oT yeseoauns oaben soon 
bas vx kextat ja pidtsakpi a S¥se ows dq Bs abvsineds B woistbaoo 
qeuds3H 3 siode in Yiaeiens og pa aah Set ove heton ite = am 
> egy? aay 15 cabada adiaaupe aie cel wilt .aodsoex g a 
vet eosses ; 
oes (Irae oye ¢(2%0 taS) ae uA pablo (doth, eo) fa : 


veied Jasiogaep Bap ftw Hiett vahneenet RS, woattoh ad? 
lea 


ae ouslg 16 beqoler 


x 


eee tape 


* ia 
ado 263 abizeaig.yisiadsta, ott sina 4¥2) 


* Hinson | ae oud ona 


98 


model. For a 3-dimensional lattice this reduces to: 


Me RE he 2 Lee, ie 


CB 


This criterion applies for an infinite lattice, and so 
it is not strictly correct when lattice truncation 
conditions are included. In effect, this limits the number 
of time steps for which the program can be run before the 
presence of the instabilities becomes apparent. FOr 
convenience, a choice of t=%/2c was made in the progran, 
and no instabilities were encountered for up to three input 
cycles. Approximately 1.5 cycles are sufficient to establish 
steady state conditions, as was verified by printouts at 1,2 


and 3 cycles. 


(5.9) Computer Program. 


A computer program based on the solution of equations 
5.9 to 5.16 was obtained from Taflove, and was re-written to 
allow it to be employed for the present problem. A _ listing 
of the new program is given in appendix B. Input data to the 
program includes a listing of the dielectric properties 
encountered at each point in a lattice of points containing 


the scatterer. The means by which this is done is detailed 


in Chapter 6. 


Before running the program with kidney data, it was run 
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for a dielectric sphere for which an alternate solution, 
using the summed-series technique of Stratton (32), was 
available. Figure 5.2 details the lattice into which 1/4 of 
a sphere was incorporated, and the parameters for which the 
program was run were €" = 4.0, 9 = 0, magnitude of incident 
Ez = 1000, f = 2.5 Ghz, 3 = \w20 = 0.3 cm, §t = 5 psec and 
envelope observation was for 460 < n < 500, where n is the 
number of time steps. For a three dimensional program, there 
is additionai flexibility in specifying the dielectric 
parameters in the free space region outside the scatterer. 
The finite difference equation for Ez (equation 5.10d) 
requires that 0(1)=0, where o@(1) is the conductivity of the 
surrounding medium. This assures that there is no 
attenuation of the incident wave. For the finite difference 
equations of E, and Ey (equations 5.10b and 5.10c), however, 
a small value of (1) can be assumed for the x and y 
directions without affecting the propagation of the incident 
wave, or the diffracted wave within the scatterer. This 
assumption results in the attenuation of the E,and Ey 
components of the exterior diffracted wave, and thus reduces 
z-directed wave reflections at the lattice truncations. This 
further aid to stability was implemented by setting ©(1)=0.1 
mho/m for eguation 5.10d and 5.10c, and setting @(1)=0 for 
equation 5.10d. The technigue could also. be applied by 
creating an H-field loss to be used in equations 5.9a, 5.9b 
and 5.9c. Judicious use of this type of anisotropy provides 


additional control over program stability. 
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Figure 5.3 compares the results of the program with 
those obtained by the exact summed series technique. It is 
noted that the computed solution locates the positons of the 
peaks and nulls of the electric field with a maximum error 
of +8, or about +3% of the diameter of the scatterer. The 
magnitude of each peak is determined with a maximum error of 
+10%. This leads to the conclusion that the accuracy shouid 
be sufficient to allow the useful solution of the fields 


within most biological scatterers. 


(5-11) Calculation of Absorbed Power. 

Since the program generates a large amount of data, 
modifications to provide results printed out in the same 
form as the lattice data were made to facilitate point-by- 
point comparison. Of primary importance in a microwave 
heating study is the amount of power deposited at each point 
in the scatterer. Because of uncertainties in the dielectric 
parameters, and to allow a compact printout, it was felt 
that a heating index in the range 0-9 would provide adequate 
information at this stage. The power absorbed in a 


dielectric medium is given by: 
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electric field, E of 1000 v/m. Normalising the 


inc @ 
dissipated power to a unit incident electric field, the 
following is obtained: 


\e1 \% 
N 2 | lEine] fra ea 


The whole lattice is now scanned for the maximun value, 
Pumas, Cf Pye and a heating index such that PRinawen en? fés 


defined as follows: 
? 
HI = N | tr eeeee SQy 


This value is printed out in the same form as the dielectric 
media lattice. The program also prints out the value of 


P so that the power dissipated at any point can be 


Niyqon 
calculated from: 
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where P;,. is the incident power density in W/m2 and is 


the impedance of free space. 
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(6.1) Modelling of Heating Processes. 


Numerical modelling of microwave heating can be divided 
into two problems. The first of these is calculating the 
fields within the body at each point, and the second is 
solving the heat diffusion equations for temperature 
distribution based on the energy absorbed at each point. The 
solution to the first problem depends on the dielectric 
properties of the body at each space point. These properties 
depend on the temperature at these points. For the second 
problem, a knowledge of the thermal properties at each point 
of the body is necessary. In a Situation where the range of 
temperatures is small, temperature invariance of the 
dielectric and thermal properties can be assumed, and a 
final steady state temperature distribution can be 
calculated by solving both problems. However, when a wide 
cCange of temperatures is involved, changes in the dielectric 
and thermal properties have to be considered. The problem is 


further complicated by changes of state. 


In the present study, microwave heating of a canine 
kidney is examined, and temperatures between -20° and +20°C 
are considered. Large changes in dielectric properties as 
well as a change of state are encountered in this range. One 


possible solution would be to divide this range into a 
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number of steps of, say, 10°C, all properties being 
considered invariant throughout the step. Considering a 
model with a homogeneous temperature distribution at the 
lowest temperature, heating could continue until the 
temperature of some point of the model had been raised by 
10°C. Based on the new temperature distribution, a second 
model of the distribution of dielectric properties could be 
derived and the process repeated. At the change of state of 
any point, a negative heating term could be introduced to 


account for the latent heat of fusion. 


Two problems arise with this approach. The first is 
that the second and subsequent models would be different for 
each new set of parameters being tested. A system where the 
new model is automatically derived could no doubt be 
devised, but the processing time, which is already quite 
long just for solving the field distribution problem, would 
become quite unreasonable. In view of these considerations, 
it is decided to limit the calculations to finding the 
heating function for various homogeneously distributed 
temperatures. Notwithstanding these limitations, several 


valuable coaclusions arise, as described in section 6.4. 


(6.2) Reasons for Choice of the Canine Kidney. 


There are several reasons why it was decided to model 


heating of the canine kidney. These are: 
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(1) Because kidneys are one of the more commonly 
transplanted organs there is considerable interest in their 
long-term storage. This has led to substantial research on 
the preservation of this particular organ, much of the 
reported work being on canine kidneys. In particular, the 
work of Rajotte (11) on the microuee rewarming of canine 
kidneys provides a valuable source of data for comparison 
with the model. 

(2) From a modelling point of view, the physical size of the 
canine kidney is attractive because a fairly detailed model 
is possible without excessive computer storage requirements. 
(see section 6.3). 

(3) Unlike some organs, the kidney contains several 
different tissue types, and this better serves to 
demonstrate the flexibility of the model chosen. In addition 
reasonable estimates of the dielectric properties of the 
various tissue types can be made from data available in the 


literature. 
(6.3) Modelling the Kidney. 


The first step in producing a reasonable model of the 
kidney was to establish appropriate planes of symmetry. 
Figure 6.1 shows dorsal and transverse views of a canine 
kidney, with typical dimensions of the organ. There are two 
planes of Symmetry, these being the dorsal and the 


transverse planes. By considering radiation along the axis 


Ls 


fo Apzneae2 Lei sansedve, a8 ane | 
eit fo Jc  yaseae ‘ 7 | 
ad¢ ,reluadsaeq ab «ayacbiy onan fron tssoury 
‘niaw 2h Ghterswen evewotole og a9 (IT) estobes R0krov 
foniisyaoo aot a7éb 20 eozioe oldsules £- eobivorg eyeshsl | . 
ee Lobos sid dain 
add to a¥dé dedteydg ode .usiv Yo 21#i0q Catlicboa 6 wort (S$) 7 : 
lajoa baltsdeb _ixist » sapezed ovldasatts ei -yembid wakaeo” 
.2nece hupet spH10ta sesuqghes av ineeaae suodtiv oldtesoq af 
-(£-@ soitoes ee) 
(79048 4iketag> yobet UAT \ensPBO sRoSe eALLau ey _ 
t wvavage qeot ai esit be ,80Qy9 ouseid -tnez920bb = 
wobtibie 0] saeaods Iebow oft to ytRigdexel2 edd saezden0med : 
<i% 0 eetztegosg Dlaroelesb, 942 Io. aadouitac eidsqcasea 
14) ab obiebtbwe £t6h aca% oben od an e601? eummhemonteee 
~orotazeehl _ 
ss we 
| — . = oer 
wdy to Isho@ @idtnaoenes 6 engonhowg ui qese gazti sdBy 8 
.vitounja 26 eSnalqg eseisqovggs wWedidayes of aby youbit 
giiaey 4 Io awety szpevadrs! OE inagob avons f 8)| ooeeet 
ows ota sted” .d6ea0 od¢ lo acofsiomkh Eeoigys dttw .yeabis 
edt bre labaodb, 243 Qnodun sade sideways to ponsiq 


nize od? Gadhia poitethn vir robiamop ga «Benntg a 


7 


— 


7 | / - 7} 


106 


EM FIELD ~___ 


| 
| 
| 
| 
| 


ele i 
DORSAL SECTION TRANSVERSE 


SE CII 


FIGURE.6.1..DORSAL AND. TRANSVERSE. VIEWS, 
WIT. DIMENSIONS, OF THE.GANINE 
KIDNEY 


SOURCE PLANE 


WIC Syhqms FAOE 


TRANSVERSE 
FACE 


20 


Soares ed 
FIGURE 62 ARRANGEMGE. OF IRE AgIICe 
AROUND AT RE Wiehe y. 


: vere ee Be 
7 i va Pp uae 


pra venaet 
Wot we 


30AA en, 


reer, SP 
- 7 ; 
P ~ 


107 


indicated on the diagram it is possible to divide the kidney 
into 4 symmetrical pieces, thus effecting a four times 
Saving in computer storage. The internal structure of the 
kidney is also approximately symmetrical about the above 
planes and so this configuration was the one selected. An 
additional advantage of this choice of symmetry is that 
radiation from above constitutes a different case from that 
of radiation from below, and so the same model can easily be 


used to examine radiation from the two directions. 


Having established the planes of symmetry, the next 
step is to superimpose a lattice of points on the scatterer 
(kidney). Experience with the test program for a spherical 
scatterer, as described in Chapter 5, suggested that 
approximately 14000 lattice points would be reasonable fron 
both CPU time and storage space considerations. Typical CPU 
time for this number of points was about 4.5 minutes with 
about 800 Kbytes of storage being used. Although 3 Mbytes of 
storage was available, the cost of each run of the progran 
was very sensitive to this factor, and it was decided to 
limit the number of lattice points at this stage. According 
to the stability criteria discussed in Chapter 5, the 
spacing of the lattice points must satisfy the following 


relationship if algorithm stability is to be assured: 
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where A, is the free Space wavelength, Amis the minimun 
wavelength in the scatterer and €',, is the maximum 
dielectric constant encountered within the scatterer. 
Ideally, the model should hold for the highest dielectric 
constant which would be encountered, in this case, probably 


that of water (#80). This gives a % (at 2.5 Ghz) of: 


=e ee Oe OG mm eke Gi 


Considering the dimensions of figures 6.1 and 6.2, the 


following number of points results (all dimensions in mm): 


N= -2 x 3° x 3 % 50,000 pts. ----:: 6-3 


This would require considerably more funds than were 
available, and so the maximum value used for €* will have to 
be restricted. In practice, this means that the model will 
have to be restricted to lower temperatures, where the value 
of €' is lower. Nevertheless, valuable information can still 
be obtained. A consideration of the required lengths of the 
various axes as dictated by the dimensions of the scatterer 
leads to the following for the number of points along each 

axis: 
nes 2A IG4 2 extra lattice spaces are 
allowed at the edges of the 


1S 
§ 
Ny = = #4 ....6S scatterer. 
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5 extra spaces are allowed at 


S the end for wave launching. 


This allows calculation of the total number of points, 
and since it was decided to limit these to approximately 


14000, the lattice spacing can be found: 
N= Noe Ny. Ng = (1243) ( #243)( 846) 2 I4 000 --- 67 


Setting k =1/6, and dividing through by 1000: 

Raita Sabet og IMERem SAO Ve 725 G'S 
ice. k 8 and $%1.25 an. 
This allows a maximum €* of 23 at 2.5 Ghz, or values of €! 
up to 80 to a frequency of approximately 1.3 Ghz. Assuming 
that @= 1.25mm, the following distribution of lattice 


points is obtained: 


N, = 15; Ny eo ees Nee ie N= 13770 


For programming convenience, an extra point was alotted 
to each of Nx and Nz for a total of 15232 points. Although 
somewhat higher than 14000, this figure was acceptable. 
Figure 6.2 shows how the lattice is arranged around the 
scatterer. The components of the impinging plane wave which 


is generated at plane J = 3 are also shown. 


Having established the arrangement of the kidney within 


- 


,2onteq to andes Less ei? to 
ylotnatsorggs of oeahs Phatl wo? bobhanmaes $4. eoatesishaa - 
sbunod, od aso paioage en4sdak od? ,O008T ie | 


‘a cane x (e uh \(e+% )(e4+3) 2 gu yh-anan | ¥ 


-0007 yd dpootds patbivkb hap ,ayr= a eatszee 
Soe + OF OFF -~ a wT + "dare «4 *A: Por 
ne 28a) TO bow OS A iawat 
‘3 to geutevy ao .2d0 2.9 a8 ES &o 'S guwizon & avolls aiat 
pnitauzed .540 Cet yiepsatxotegs yo yoredpett 4 of 08 oF Qu 
svldssl 26 soldadiatets pakvolget od .aa@Sst =F pads >Hi 


OMets BM  fTSom gk \ et = % scr, = 4s ay a0 
owe 
=) ae 

bedtole. ase zee Syne os. .Soeelcvuvuss pPHleasipoTg 3OL . a” | 
ApuodelaA .adtaiog, SESE Re Leto & 33 gh Dae nt io. sone) +081, 
«9idsiyeons die onupi2? eft3 .000RT sads 2adp di Sear) 4 
oi) hovers bepastys 22 soissnt eto gon avoda S.8 emop.a | 
torsy ovse sakg paipdicat odd to 2tnegeqe@oo od? cae 7 
-heode Seles ome & * G paalg 24 ae ak 


i402 yonbix ody Fo duwnegnsrse Odd Pi: MRD oth a 


‘a 
. ary 7 . 


the lattice, it is now necessary to establish the internal 
details. To minimise the number of pages of data printout, 
it was decided to construct 15 planes (K = 1 to K = 15) of 
27x34 (IxJ) points. This necessitated visualising the kidney 
structure at each plane. To aid in this procedure, figures 
6.3, and 6.4 were obtained from standard texts of canine 
anatomy (58,59). Using these drawings, and designating the 
different tissue types as indicated below, the set of planes 
of figure 6.5(1) to 6.5(15) were drawn. The bracketed number 


refers to the plane number. The tissue types are: 


1 - medium external to scatterer, in this case, air. 
2 - cortex 

3 - medulla. 

4 - renal siaus (fat). 


5 - renal pelvis (connected to ureter). 


Although figures 6.3 and 6.4 do not make the disposition of 
the renal sinus clear, Evans (58) describes ant as 
"Surrounding the renal pelvis", and this is how it is 
considered here. Many details are too small to appear on the 
1.25mm lattice. This includes the vascular structure, 
comprising, among other things, the artery and vein leading 
to the kidney. Although these are in the same size range as 
the ureter, the latter becomes the relatively large renal 
pelvis within the kidney, while the vascular structure 


becomes progressively smaller as it evolves into the 
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FIGURE 6.5(1-15) 
The 15 planes of the dielectric media lattice. The media 
are: 
(1) Air 
(2) Cortex 
(3) Medulla 
(4) Renal Sinus 


(5) Renal Pelvis 


Note: The labelling of the I planes is somewhat unusual 
and is as follows: 
On figures V2 a awiemicc cise 1a IB Ue oc cea - eke ee 
Oe le 2 56°77 
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capillary system. Because of the importance of the vascular 
structure in the perfusion of kidneys during the freezing 
process, it is hoped that some means will be found to 
incorporate this structure into the model. The limitation 
imposed by the exclusion of this structure is that 
alternative means must be found to account for the presence 
of the perfusate within the other parts of the kidney. This 


is done by appropriate choice of dielectric properties. 


SSS See Se SS Se SSS qe a Se ee SSS ee Se Se eS ee 


In order to make use of the model, it is necessary to 
assign dielectric values to the various media. The only 
known source of data for the properties of DMSO-perfused 
kidney are the data of Foster et al (22), and this data 
applies only to media 2 and 3 at 1 Ghz. In addition, the 
degree of perfusion of the various media is not definitely 
established. Therefore, estimates based on the existing 
literature data and assumptions regarding the degree of 
perfusion of the various media must be made in order to 
assign suitable values for the model. At this stage, the 
model can be used to make predictions in two important 
areas. These are (a) in establishing the amount of energy 
absorbed by the kidney at low temperatures, thereby allowing 
some estimates of possible heating rates, and (b) the effect 


of different parameters on heating uniformity. 
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It was decided to examine the two above areas in the 
light of four different parameters. These are: 
(1) Illumination frequency: The two ISM bands primarily used 
for microwave heating are at 915 Mhz and at 2450 Mhz. 
Illumination frequencies of 1000 and 2500 Mhz were used in 
the model because much of the existing data pertains to 
these frequencies. 
(2) Cryoprotective Solutions: Most reported studies of the 
cryopreservation of kidney tissue have used 10% DMSO. For 
purposes of comparison, it is interesting to examine the 
effect of a 1% DMSO perfusate as well. 
(3) Temperature of model: For reasons outlined earlier, it 
is difficult to model a continuous heating process, and so 
it is necessary to choose a particular temperature at which 
the model is to be tested. At 2500 Mhz, the maxinun 
temperature at which the model can be used is limited by the 
maximum dielectric constant encountered. To give some idea 
of how heating patterns would change with temperature, it is 
interesting to examine the model at different temperatures. 
The choice of -20°C and -10°C was most reasonable at 2500 
Mhz. At 1000 Mhz, there is no temperature limitation, and 
so, in addition to the choices of -20°C and -10°9C, a 
temperature of +20°C was also selected. The use of the model 
at temperatures below -20°C was precluded by lack of data. 
(4) The kidney state: 4 different kidney states are defined 
in this case. Foster's measurements (22) and the relevant 


measurements based on frozen foods and other data, as 
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presented in tables 6.1, 6.2 and 6.3 were used to obtain 
appropriate estimates of the dielectric properties of the 
Various media. Table 6.4 shows the properties used in the 
model as estimated on the basis of the discussion below. The 
four kidney states are: 

(a) No perfusion: This state was chosen for purposes of 
comparison, and is also useful for examining the effect of 
the perfusate on the heating rate. 

(b) Selective perfusion: Most perfusion studies suggest that 
the cortex and medulla, which have an extensive vascular 
system, are well perfused. The pelvis, which contains urine, 
would certainly undergo some perfusion by virtue of the 
penetrative properties of DMSO. The renal sinus, on the 
other hand, is comprised of fat, a medium in which perfusion 
appears to be minimal. For the present work, a well perfused 
cortex and medulla, a partially perfused pelvis and an 
unperfused sinus were selected. 

(c) Pelvis completely filled with perfusate: This is an 
unrealistic assumption from a physiological point of view, 
but because of the uncertainties of the dielectric 
properties of urine, it was decided to assume that the 
pelvis was completely filled with the perfusate, and would 
therefore be more lossy than assumed in (b). Because some of 
the tubular structure of the pelvis appears in the model, 
this condition would also serve to give an indication of the 
behaviour of the vascular system, which normally also 


contains perfusate. As previously explained, the vascular 
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structure does not appear on the model. 


If all possible combinations of the above parameters 
were to be modelled, the program would have to be run 36 
times, and even if appropriate funds were available, 
meaningful analysis of the data would be an enormous’ task. 
Therefore, a datum set of parameters was selected from the 
above possibilities, and selected changes from the datum 
were made in order to examine the effect of the other 
parameters. The following set of conditions were chosen as 
the datum because they represent the situation most likely 
to be encountered in the real world. The conditions are: 
(a) Illumination at 2500 Mhz. 
(b) 10% DMSO perfusate. Realistically the perfusate would 
also contain some ionic components, such as NaCl, but 
reference to figures 6.17 and 6.18 suggest that these have 
little effect on the dielectric properties. 
(c) A temperature of -20°C. This would be the case with the 
lowest loss and would be most representative of the low 
temperature heating rate. 
(d) The partially perfused model. This is most consistent 


with the measurements reported in the literature. 


Table 6.4 shows the 8 different sets of parameters to 
which the model was applied. The dielectric properties 


assigned to the media in the table were calculated on the 
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basis of the following assumptions: 


(1) The frequency dependence of the cortex and medulla was 


anaes all... ha la= tt a i oo ae ik Se he ee 
| Material |€"]| Freq. | -20°C | -10°C ] +#20° i Ref. | 
| ini e Shhze E! | { | | 
-——___——_+—_ +--+ ++“ 
| {E*] 900 | 4.8 * | 6.2 * 458.0 * |] 46 | 
{ Beef 17 | { 0.046 | 0.103 | 1.226 | { 
| (Eh) eo. 27008 [R428 *£ehio6s 2'42(47.7 e939 | 
| io | ieO. 079s (SOs 210 9 | hoe 63R i { 
Ss a a: Sa Sea Sen meee | 
| $e"; 1000 | 2.0 * | 3.3 * J60.0** | | 
{Unperfusedj{c | [--05002-j~ 0..0084|—1..250-«1).-22+-| 
ji Kidney jé€'j 2500 j 2.0 e373 {50.0 | | 
f lo | {| 0.003 { 0.017 | 1.900 | | 
a 4 
| 1% 4E€"} 1000 jf 2.1 [esos {60.0 | { 
{ Perfused |o | [702 00SS 10. OTR et. 2209} { 
[-—Kidney- 6 [-E%|--—-2500——{-2..1 { 3.5 {50.0 | | 
| lao } | 0.004 | 0.023 { 1.86 | | 
[ge ee a a pee SS a ee | 
| 10% ie"{ 1000 {4 3.0 * { 5.0 * {60.0 * | | 
| Perfused |o | [ 0. G099[(* 0.035%] 1.000 | 2271 
{ Kidney jé€*j 2500 | 3.0 fe5.0 {50.0 | | 
| lo | jn0:0:0 15 e¢[+02:078-—|—12520—| { 


Notes: (1) Vaiues marked * are from the literature, 
and all others are calculated. 
(2) Values marked ** are similar to Schwan's 
(19) 


TABLE 6.1 Dielectric Data for Beef and Kidney 


assumed to be the same as that of raw beef, for which data 
at both 1000 Mhz and 2500 Mhz were available, as shown in 
table 6.1. In addition, the frequency and temperature 
dependence of perfused and unperfused kidney tissue was 
assumed to be the same for the small temperature ranges 
considered. The effect of the perfusate on the tissue 


properties was assumed to be linear with concentration. The 
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values marked with an asterisk in table 6.1 are the values 
measured by the investigators indicated in the reference 
columna, and the other values were derived from these. 
Referring to table 6.1, note the large difference between 
the values quoted in the literature for the conductivity of 
beef and kidney at low temperatures. This is surprising in 


view of the similarity of their properties at +20°C. 


eee = 7 1 
{ Material)€*"| Freq. (—209Cee T= 10°C Fe] 420°C ee P43 7 8C Si ker. 
| loi Mhz. | | | | l | 
eee eee 
i Lard j€"]j 2500 | [228 ee [2556480] 2268xe5| | 
| (Pig fat) |o | | {0.006 40.013 [0.018 439 | 
SS pe 
| Fatty J€'"| 1000 25.5 225 PA eS) | i | 
| Tissue |o] {0.008  |0.016 40.032 jf | | 
Ey ee a ne | 
j Fatty j€*| 3000 Ze PAS) [2Ze.5 | H { 
| Tissue |o|] 10.006 40.012 40.024 | { | 


(bee ed el Ed lh 


Notes; (1) Values marked * are measured or literature 
values, the others are calculated. 


TABLE 6.2 Dielectric Data for Fat 


(2) Literature data for fat varies considerably, and there 
is not much data available for low temperatures. However, 
assuming that the temperature and frequency dependence of 
different fats is the same, the vaiues in table 6.2 were 
estimated on the assumption that they would behave in the 


same way as lard, and using the literature values indicated 


with an asterisk. 
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Le 2 eS I Ce Seen a aay Kare mey nay [aan eae Eee 
{ Material J€"| Freq. 4-20°C 4-10°C 4#209C |Ref.| 
l iPtlisathz. *| | { | i 

rae ee 
j FES fe 1000% [= 3e2 {| 3.2 * {79.0 * | | 
J Water lo | 0.001 { 0.001 { 0.280 | 43 | 
| fe! Pee2500° efy3s2 { 3.2 * {78.0 * | | 
| lo | {| 0.001 { 0.001 | 1.530 | | 
a Pr a Se Se Sa cae 
| 1€*] 1000 | 2.9 jaesou3 {75.0 | | 
| 1% | {| 0.040 | 0.059 | 0.290 | { 
{ DMSO .{€"fgpe2700 [ 2i98*ee(53289* 173.5* och! 
j Ic] | 0.040 4 0.059 |] 1.558 | 6 | 
Hh th eH A 
| 1€*f 1000 | 7.5 {18.5 174.0 | | 
| 10% Io | {| 0.652 | 1.830 { 0.380 | 
i DMSO JE" ] Z2700° | 7725 * 11855 173. 0*OgisCh. | 
| ica | { 0265291105830 |.2.05280) 6m 
(dd a a te ee 


Notes: (1) Values marked * are measured or literature 
values, the others are calculated. 


TABLE 6.3 Dielectric Data for Water and DMSO 


(3) The values of table 6.3 are derived on the assumption 
that the behaviour of DMSO solutions for low concentrations 
is similar to that of water, both as regards frequency and 
temperature. This is justified by the results obtained in 
Chapter 4 Assuming 10% of the volume of the urine in the 
pelvis is displaced by perfusate, and assuming linear 
mixing, the values for the pelvis are calculated and shown 
in table 6.4. It is assumed that the properties of water are 
Similar to those of urine for temperatures of -10°C and 
-20°C. Again, this is justified by the results of figures 
6.17 and 6.18 where an ionic component has little effect on 
the low temperature dielectric constant. In any case, most 


of the losses are due to the DMSO. 
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NREL TE RSG) eR aaa, PAN RAT ae poems SRR MRS Sa rs gee gE rE | 
{Cond.{ Freq. | Per- {| Temp.| Kidney | Dielectric |Remarks| 


1 # | Mhz {|{fusatej °C | State |{ Properties | | 
Hh tH He TI 
| | | | J1- 1.0,0.005] | 
| | { 10% | { Partialj2- 3.0,0.015j| Datum | 
1 1 | 2500 § DMSO | -20 | per- 1|3- 3.0,0.015] model | 
| J | | { fusion |4- 2.5,0.006| | 
| { | | | {5- 3.6,0.066| | 
-—_+—_——__+—___—_+-——- +11 
| | | | | {1- 1.0,0.005] | 
j | 1 10% j4| | Partial|2- 3.0,0.009] Freq. | 
He PA | 1000 | DMSO {| -20 { per- {3- 3.0,0.009{ effect] 
I | { | | fusion |4- 2.5,0.008] | 
| i | | | [5- 3.6,0.066] | 
r TE eT | 
| | | | { }1- 1.0,0.002| | 
| { i 1% | ] Partial|2- 2.1,0.004] Per- | 
1 3 4 2500 | DMSO y -20 | per- {3- 2.1,0.004| fusate] 
| { i | {| fusion {4- 2.5,0.006]| effect | 
| | | { | {5- 3.1,0.005| | 
i I ree 
| | | | | [1- 1.0,0.005] I 
| | 4 10% | { Partial|2- 5.0,0.074j Temp. | 
i 4 qj 2500 j; DMSO jy -10 | per- {3- 5.0,0.074j effect| 
| i { | | fusion {4- 2.5,0.012] | 
i | | | | {5- 4.7,0.180} i 
| SS eo aa IS eee oon fees ae pe gee aS pe eS | 
| | | | | [1- 1.0,0.001| | 
{ { 4 10% | | No }2- 2.0,0.003] Per- | 
{ 5S 4 2500 | DMSO {| -20 { per- {3- 2.0,0.003{ fusion| 
] j { | | fusion {4- 2.5,0.006] effect| 
| | | | [5- 3.2,0.001| 
Se a ee ee 
| | j | | Total j{1- 1.0,0.005] Effect] 
{ { {10% | { pelvic |2- 3.0,0.015] of | 
}°°6 °j°°2500 7} DMSO | -20 } “per- [3- 3.0,0.015| lossy | 
| | | | {| fusion |4- 2.5,0.006] pelvis| 
r |5- 7.5,0.650| 


-—_—__t—__—__t—___-—_+————_t+__—_1——q—_I1 
J}1- 1.0,0.100]{ Temp. | 


l | | l | 

| { 1 | { Partial|2-60.0,1.000] and | 
i 7 | +1000 | DMSO | +20 | per- {3-60.0,1.000{ freq.| 
| | | | | fusion |{4- 2.5,0.032| effect | 
| | { | | {5-78.0,0.290] | 


| 1.0,0.005] { 
{ j Partialj2- 3.0,0.015jOrient- | 
Oj; -10 j per- j3- 3.0,0.015{ ation | 
| { fusion |4- 2.5,0.006] effect] 
| | {5- 3.6,0.066]| | 
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(4) fhe particular choice of the assumed conductivity of air 
is to keep losses in the air to a minimum in comparison with 
those of the kidney, and yet to ensure algorithm stability. 


The values chosen appear in table 6.4. 


eS es Ss awe Se 


The eight sets of conditions for which the model was 
run generated 120 pages of data printout and it is obviously 
impractical to present all of these. It was therefore 
decided to consider each set of conditions separately and to 
present oniy those planes which showed important results. 
Where appropriate, comparisons between the different sets of 
conditions are drawn. General conclusions and comparisons 
with the experimental data of Rajotte (11) are presented in 


a concluding section. 


As explained earlier, the printout is arranged in 
exactly the same manner as the media lattice input. A digit 
in the range 0-9 is printed out, and this represents the 
heating index at each lattice point. The scaling factor to 
obtain the actual value of power absorbed in W/m3 at each 
point is indicated in each case. This value is for unit 
incident electric field and can be used to obtain the power 
absorbed for any incident power density by means of equation 
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For all the cases where the incident frequency was 2500 
Mhz, the program of appendix B was run for 480 steps with 
printout at 400 and 480 steps to enable a stability check to 
be made. No instabilities appeared for any case and only the 
results obtained after 480 steps are presented here. As a 
result, the program was run 600 steps in the 1000 Mhz case 


with only one printout (at 600 steps). 


For each case, the dielectric properties are presented 
in table 6.4. The conditions under which the dielectric 
properties were obtained are also shown there and will only 
be repeated where necessary. In labelling the printouts, one 
figure number at the beginning of each set will be used to 
represent all the planes within that set, with a bracketed 
addition to refer to a particular plane within the set, as 
was done earlier. For instance, the 10th plane in figure 6.6 


will be referred to as figure 6.6(10). 


(6.6) Condition 1. The Datum Model. 


(Figure 6.6 (8-11, 13, 15) 


This figure is the printout obtained with the datun 
model. Plane 8 shows that the dissipation in the medulla and 
cortex is minimal, and planes 3 to 7, which are not shown, 
support this observation. The maximum index encountered in 
any of these planes is 3. In plane 9, an interesting effect 


becomes evident. Here it is clear that the tubular structure 
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FIGURE 6.6(8-11,13,15) Datum Model 


Set #1: Max. dissipated power = 0.0273 W/m3. 
Max. dissipation at plane K = 9. 
Frequency = 2500 Mhz. 

Perfusate = 10% DMSO. 
Temperature = -20°C. 


Partial Perfusion. 
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NORMALISED HEATING INDEX AT TIME STEP 480 
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of the ureter is dissipating strongly in contrast to the 


surrounding fat layer in which no power is being absorbed. 


damage__due to heat stress and mechanical stress_caused_by 
differential heating. Referring to planes 10,11,13 and 15, 
it is clear that this effect becomes gradually less 
pronounced as the tubular structure becomes larger. However, 
in none of the planes does the fat absorb significant power, 


and the fat-pelvis_ interface is___always___subject___to 


considerable differential heating. 


Further reference to planes 10-15 indicates that most 
of the energy is absorbed by the perfused urine in the 
pelvis, in spite of the supposedly good perfusion in the 
cortex and medulla. The best heating in the cortex and 
medulla occurs on the side away from the impinging wave. The 
effect of the small loss introduced into the air is seen to 


be negligible in this case. 


(6.7) Condition 2. Effect of Frequency. 


(Figure 6.7 (8-11, 13,15) ) 


The first effect observed in changing the frequency 
from 2500 Mhz to 1000 Mhz is a drop in the maximum absorbed 
power from 0.0273 to 0.0117 W/m. Again, planes 3-8 (only 8 
is shown) in the off axis region of the cortex and medulla 


show little dissipation. Some loss in the air is evident at 
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Set #2: Max. dissipated power = 0.017 W/m3. 
Max. dissipation at plane K = 11. 


Freguency = 1000 Mhz. 


Perfusate 10% DMSO. 
Temperature = -20°C. 


Partial Perfusion. 
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this frequency, but this is due to the artificial loss 
introduced for stability and only those regions within the 


kidney should be considered here. 


Planes 9,10 and 11 indicate a similar tubule heating 
effect as was observed for condition 1. Planes 11,13 and 15 
show that the pelvic region is again strongly heated with 
even less heating in the medulla and cortex. At this 
temperature, the higher frequency (2500 Mhz) appears to give 
Slightly better uniformity and approximately doubles the 


maximum dissipated power. 


(Figure 6.8(5,8,9,10,12,15) ) 


Again, it is immediately evident that the maximun 
dissipated power is reduced, in this case by a factor of 
more than 6. This bears out the prediction made in Chapter 4 
with regard to the loss factors of 1% and 10% DMSO. With 
regard to the heating in the medulila and cortex in this 
case, it is apparent that a much higher proportion of the 
energy is absorbed in this region. Planes 3-8 all show 
reasonably high absorption, with a maximum heating index of 
9 in plane 5. Only planes 5 and 8 of this group are shown 
here. Because it was necessary to set the air loss 
relatively high in this case, it is clear that some energy 


is being dissipated in the air, and so, as before, only the 
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FIGURE 6.8(5,8-10,12,15) Perfusate effect 


Set #3: Max. dissipated power = 0.0033 W/m3. 
Max. dissipation at plane K = 5. 
Frequency = 2500 Mhz. 

Perfusate = 1% DMSO. 
Temperature = -20°C. 


Partial Perfusion. 
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losses within the kidney should be considered. 


Although plane 9 does indicate some differential 
heating, it is the fat layer which is more dissipative in 
this case. Further, although the highest dissipation is 
still in the pelvic region, dissipation in the cortex and 
medulla, as indicated by planes 10,12 and 15, is a much 


higher percentage of the maximum. A reduction in the DMSO 


uniformity, although it should be remembered that_this_is_at 
one __ particular temperature only. The relatively small 
variations in conductivity in this case suggest that more 
uniform perfusion at any concentration, thereby resulting in 


more homogeneous dielectric properties, would certainly 


enhance heating uniformity. 


(6.9) Condition 4. Temperature Effect. 


(Figure 6.9(5,8-11, 13, 15)) 


In increasing the kidney temperature from -20°C to 
-10°C, an increase in maximum dissipated power of more than 
three times from 0.0273 to 0.0920 W/m3 is observed. Planes 
3-8, of which 5 and 8 are representative, suggest that 
somewhat improved heating in the cortex and medulla has been 
obtained. Again, the tubular effect at plane 9 is evident, 
but planes 10-15 suggest that dissipation in the pelvic 


region is relatively lower and overall uniformity is better. 
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FIGURE 6.9(5,8-11,13,15) Temperature effect 


Set #4; Max. dissipated power = 0.0920 W/m3. 
Max. dissipation at plane K = 10. 


Freguency = 2500 Mhz. 


Perfusate = 10% DMSO. 
Temperature = -10°C. 


Partial Perfusion. 
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Maximum dissipation occurs in a tubule in plane 10. 


Although conclusions are difficult because heating 
would not be from a uniform temperature in the real case, it 
seems that a slight improvement in uniformity is obtained as 
the temperature increases. However, dissipation in the fat 


layer is still minimal. 


(6.10) Condition 5. The Effects of No Perfusion. 


a Se Se a Se ws Se ee Se es Se a Se Se Se eee 


(Figure 6.10(5,8,9,11,13,15) ) 


As Might be anticipated, the maximum dissipation is 
considerably reduced, to approximately the same level as for 
1% DMSO. Dissipation in the cortex and medulla, of which 
planes 5 and 8 are representative, is relatively good, and 
the tubular effect at plane 9 is reversed in that the 
tubules themseives absorb no power. In contrast to the 
Situation for 1% DMSO, however, it is clear from the 
representative planes 11,13 and 15 that no dissipation is 
occurring in oye pelvis. In view of the sensitivity of the 
conductivity of this region to the DMSO concentration, this 


is not surprising. 


(Figure 6.11(9-11,13,15)) 


In this case, the conductivity of the pelvis is 
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FIGURE _6-10(5,8,9,11,13,15) Perfusion effect 
Set #5: Max. dissipated power = 0.0032 W/m3. 

Max. dissipation at plane K = 15. 

Frequency = 2500 Mhz. 

Perfusate = 10% DMSO. 

Temperature = -20°C. 


No Perfusion. 
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FIGURE 6. 11(9-11,13,15) Tubule effect 
Set #6: Max. dissipated power = 0.254 Wym3. 
Max. dissipation at plane K = 9. 
Frequency = 2500 Mhz. 
Perfusate = 10% DMSO. 
Temperature = -20°C. 


Partial Perfusion with 10% DMSO in pelvis. 
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increased to a value representative of total displacement of 
the urine by the DMSO perfusate. This is interesting for two 
reasons. The first is that it further indicates the 
region, and the second is that it suggests what might happen 
to the vascular system, since this is likely to contain a 
very high percentage of perfusate, perfusion being through 


the vascular system. 


Plane 9 is the first plane to indicate any dissipation 
whatsoever, and the value here is very high, and is in the 
tubular structure. Maximum dissipation is 0.254 Wym3, a 
factor of more than 9 times that of the datum model. At this 
plane, there is no absorbtion anywhere else. Planes 10,11,13 
and 15 are representative of the rest of the kidney, and 
show that all the power is absorbed in the pelvic region. 
From these results, it is clear that a_high content of 


perfusate in smali_ tubular _ structures is liable __to__cause 


Sa a ——. me Sas ee Se 
ee ee ee ee a wea SS 


potentially damaging results. This suggests the possibility 


of a final perfusion step with an appropriate liquid chosen 
to minimise this effect. 


(6.12) Condition 7. 1000 Mhz Heating at _20°C- 


(Figure 6.12(5,8,9,11,13,15)) 


The first interesting point to emerge from this run is 
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Set #7: Max. dissipated power = 0.115 W/m3. 


Max. dissipation at plane K = 10. 


Frequency 1000 Mhz. 
Perfusate = 10% DMSO. 
Temperature = +20°C. 


Partial Perfusion. 
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that the maximum dissipation under these conditions is only 
Slightly more than that obtained at -10°C and 2500 Mhz 


(0.1150 W/m3 as compared with 0.0920 W/m3) This supports the 


earlier observation that more energy is dissipated at 2500 


Mhz, and therefore heating rate is likely to be better at 


frequency. Planes 3-8, of which planes 5 and 8 


are representative indicate reasonable absorption in the 
off-axis planes of the cortex and medulla. Plane 9 indicates 
that there is no heating in the tubules or in the fat layer 
surrounding them. Planes 11,13 and 15 are representative of 
the axialiy located planes, and they suggest that very 
little energy is dissipated in this region, with none at all 
being dissipated in the pelvis and sinus. This is in 


contrast to the situation existing in most other models. 


(6.13) Condition 8. Reversed Radiation Direction. 

No figures are presented for this condition, but 
results showed that there was little difference in the 
behaviour of the model for this condition as compared to the 
datum model. This may not necessarily be true for higher 
loss conditions, but it does suggest that there is no 


particular sensitivity to orientation. 


(6.14) Comparison with Experimental Results. 


Rajotte (11) performed a number of microwave heating 
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experiments utilising canine kidneys. A brief summary of his 
experimental conditions and results follows: 

(a) Freezing: After rinsing in Baxter's solution, DMSO was 
added in a two step perfusion process at either 4°C or 25°C. 
In the first step, 10% DMSO in Ringer's solution was added, 
and in the second, 20% in the same diluent. Cooling was by 
liquid nitrogen cooled fluorocarbon perfused through the 
vascular system. The cooling rate was controlled by the 
differential temperature between the ureter and the cortex. 
(b) Thawing: Thawing was in a multimode cavity at 2450 Mhz. 
The kidney was immersed in a teflon container of 
fluorocarbon, which, like teflon, has very low loss at 
microwave frequencies. After thawing, internal temperature 
gradients were measured with a thermocouple probe. No 
gradients were found in the cortex, but cold spots were 
found in the medulla and pelvic region. Adequate perfusion 
resulted in temperature gradients of +10°C but inadequate 
perfusion gave much larger gradients. Urine remaining in the 
pelvis caused a large temperature gradient. When the kidney 
was cut in half, following thawing, it was clearly seen that 


areas of incomplete perfusion were still frozen. 


Direct comparison between the experimental results and 
those of the model cannot be made due to the entirely 
different illumination conditions. However, some general 
comparisons can be made. In the first place, the nodel 


confirms that gradients are to be expected if urine remains 
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in the pelvis. A comparison of conditions 1, 3 and 6 
Suggests that a lossy pelvic region would strongly dominate 
the overall dissipation pattern, and it is obviously 
advantageous to have a fairly low loss material in the 
pelvic region. This would also aid in reducing the tubular 
differential heating which was observed in most cases. It _is 
bossible that careful adjustment of the loss_in this region 
1,3 and 5 suggest that there may be an optimum DMSO 
concentration (up to 70%) for a given temperature, and 
techniques__whereby__concentration is_varied during heating 


may_also_ provide a_means of controlling uniformity 


# 


The model confirms that gradients in the cortex are 
minimal, although it suggests that heating in this region is 
hot very strong. Rajotte's (11) experimental results do not 
bear this out, but the fact that heating in this case was in 
an oven, where radiation impinges from many directions, and 
was also at temperatures up to 20°C, could explain this 
discrepancy. Evidence of the effect of temperature is 
provided by condition 7 where good heating of the cortex and 
medulla were obtained with little dissipation in the pelvic 
region. Though the frequency is different, this behaviour is 


Similar to that found by Rajotte (R8). 


Further consideration of condition 7 shows that there 


was little absorption in the pelvic region and tubules for 
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this run, in spite of the relatively high value of 
conductivity. These findings agree with Rajotte (11) whose 
experimental results indicated cold spots in the pelvic 
region after heating. The model assumed poor perfusion for 
the pelvic region and so it seems necessary to ensure better 
perfusion if adequate pelvic heating is to be obtained. This 
also supports Rajotte, who concluded that the pelvic region 


was often poorly perfused. 


The experimental measurements of temperature gradient 
were rather coarse, and only took place after the kidney had 
been thawed. It is therefore very unlikely that the 
gradients at the fat-pelvis layer would have been observed. 
Other investigators (53) have found similar gradients in 
other situations which suggests that those observed in the 


kidney model do indeed exist. 


(6.15) Concluding Remarks. 

In conclusion, the model has revealed several important 
aspects of microwave heating of perfused canine kidney. {In 
the first place, addition of DMSO as perfusate increases _ the 


maximum dissipated power, and this would certainly lead to 


increased heating rates at low temperatures. In this case, 


more energy was absorbed at 2500 Mhz than at 1000 Mhz. 


Secondly, it is clear that uniform perfusion is 
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essential if uniform heating is__to be achieved. This 


requirement is in conflict with the present methods of 
perfusion whereby almost pure perfusate is likely to reside 
in the vascular structure. This results in highly non- 
uniform heating, especially in those regions which are 
adjacent to a fatty layer. This condition is further 
aggravated by the thermal insulation properties of fat (23), 
which would tend to reduce heat conduction away from the 


strongly heated areas. 


uniformity _is_ quite dependent on the properties of _the 
pelvic region, and appropriate choice of properties in this 
region may provide a means of controlling uniformity. The 
model also suggests that maximum uniformity at a given 
temperature could be obtained by appropriate choice of DMSO 
concentration. Unfortunately, adjustment of DMSO 
concentration in the tissue is difficult, but it could also 
provide a means of uniformity control. DMSO perfusate for 
best uniformity at each temperature, and continuous control 


of concentration may also provide a means of uniformity 


control. 


Finally, the model assumed that there was no perfusion 


in the sinus, and as a result, almost all of the different 


conditions modelled indicated sharp gradients in this area. 


=e ee es ee ee ee ee ee ae we re ee 


ehiees oe 


ora. ded. :iallensdialine ak 
-adjiad 22 moddhbnos) skit aayol ihn 2 omeeths 
(en) vot 2a pebssaging uottatvaai Lone do Ua bedewnaeye 
vas @omt Yshe aoktoahaos —s6aH ojunee @¢ Bast Sinow toby 

.2a7s bossed yLguorse 


. wht san 
liaseyo tedt. ob azepisas Jed) 20t56% sasdrogns biidt Ao 
shtt ak aetaasqedg Io coicds ddd sqooqqa hae siobpes _pvivieg 
od? .ytteae®iay paiilorteep Jo aasem 6 ebivcrg yan molges 
oovie ots (obnzol ied ovetxsm, todd atespeve onls feboa 
oend to eododo etsitqotg9s Yd fsntetdo at binoe sisssieqgas? 
of ad 10 tgousaut he VY fogedy 4 totau -colts13 aeDnta 
oels filons ttl gud ,sLlrotthit 22 oneets end? at sol¢sa? GeDao0 ; 
tok etouplaed WRRd stor 7aen \Yeiasotbew Ja enson 8 “obkvoag 
lovtaes sauouakties ba» .smwtesegqess dose 3a \Jiazotiay seed 
yrimpotios Jo este 5 obiworq oaks ype sBolt529nN9DR0R'~ RO 
-Loiss0o 
wee 
aoiewhay on sev o2edt edd bomvers Lebom edt .yiisoti ©+% 
thoretedh off 10 fhe toons ,tiveez a a6 bas .euoie emt Re 
.eoas ebit af efaekbete quaie bsteolbai belipboa anotsibaoo 


useful in checking this result, and_means Oo __improve 


it is found to be lacking. 


Regarding accuracy, the results obtained with the model 
would certainly be improved if more accurate estimates of 
the dielectric properties of the various media were 
available. However, the conclusions which can be drawn at 
this stage already demonstrate the utility of the modelling 
technique. Other improvements which would expand the 


usefulness of the model are discussed in Chapter 7. 
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Chapter_7. Suggested Further Research. 
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The application of the numerical model of Chapters 5 
and 6 indicated a lack of appropriate data for use in 
numerical modelling of microwave heating in cryobiology. 
Extension of the model to considerations of heat conduction 
wili also be hampered by the lack of suitable data for the 
thermal properties of tissue. In view of the large amount of 
data required, automated technigues for making thermal and 
dielectric property measurements are therefore very 


desirable. 


The technique for measuring dielectric constant as 
employed in Chapter 4 is amenable to automation. The use of 
a microprocessor is appealing. The processor could adjust 
the frequency of the source for maximum transmission through 
the cavity, calculate the frequency shift, note the 
attenuation, and then, after applying appropriate 
corrections, calculate the dielectric constant. In 
conjunction with the temperature calibration curve, the data 
could immediately be plotted on an x-y plotter to give 
dielectric constant vS. temperature curves. Extension of the 
model to other frequencies would require an expansion of the 
dielectric measurement facility. Because most of the other 


frequencies of interest from a heating point of view are 


198 


@ paesqeia Xo febod naan he 


a1 eeu 20% atoh, etebaqoagqe” X oet we beraodbak 9 bas, 


-\eeloidoyz of gakipen dvshosotn (20 ablichos Lepineaum rc 


uoltaybaon ‘wot to doi teo4bkadsp of iehow @d3)2o -actenesat 
ad? sot egph elubsive te domh ele poliemegend od woteuiitte ” 


10 ¢apoas apisl ed? rc velv, OF usvante 29 esitzeqozg Lanredt - 


bes £saxeds pnadae (6% aompiodged> HOT SHORE © ,beaiupsa sd6b 
ytev oaeteresd o15 Sf neste nism yaueqosq oigzveietb — 


or 


as tustseo> otadpolsib patzaessa 702 supiadoas eit. 
eo «say of .nobeewotus, of siicatws ak * se2qsd. ah beyoigas 
teutis biseo tOaeeter of! -pnkiseggn af tozegn07gora ke #) _ 
A:HO7A% acres iaBnees avaixed 761 SOIR ad¢ Yo youespe7? oid 
aé+ of00 <—h3id2 Yousnpavt sd% eradupiss .\siveo. say 
ocwlsgeiqys gnayiggs setts ~ .wodt. bas cohssuassss 
al -tastanago Bixfosioib sda et¢éieviss ,anolsoezi0> 

cteb odd? yovie aobdntdkies eastezeqaed Gado dite ashes ae 


avip of nettolg y-x 40 so bettolqg od ylossibonat biso> 


edt to moteaosed .aevi09 etu¢s2eqast av tastanoo vizsoelelh 
ait to ackensqxs de siivpe2 bigow asinasypeni r9d30 of Lebou 
wsito 9nd to taou esunned “VIiLos} dovassvesem ob asaeiedb 
#xs vote Yo. takoy yuitsod & sox gaszessk to neinseppead 


7 


199 


lower than 1 Ghz, the use of the time-domain method of 
Iskander and Stuchly (35) is attractive. Using this 
technique, a single measurement suffices to give data for a 


wide range of frequencies. 


Regarding thermal properties, experience with frozen 
foods (23) has suggested that the best approach is to use 
the enthalpy of the various tissue types. This would provide 
a measure of the total energy required to raise the 
temperature of a particular tissue by a given amount. When 
used in conjunction with thermal conductivity data obtained 
under identical conditions, enthalpy data would thus allow 
calculation of the temperature distribution within a body 
subject to microwave or conductive heating. In addition, if 
emissivity data were available, the effects of conductive 
and radiative heating could also be modelled, and the nost 


effective combination could be determined. 


(7.2) Numerical Heating Model. 


The model described in Chapters 5 and 6 represents a 
first attempt at modelling microwave heating processes in 
cryobiology. In the present work, only a heating index for 
different model parameters was calculated, and general 
conclusions as to the effect of plane wave radiation were 
drawn. In an actual heating situation, however, it is the 


resultant temperature distribution after heating through a 
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range of temperatures that is important. In addition, plane 
wave heating is rarely encountered in practice, and 
application of the nodel to heating in a microwave oven 
would be desirable. Increasing the resolution of the kidney 
lattice to allow incorporation of the vascular structure 
would also be an important improvement in the model. 
Finally, extension of the model to an examination of the 
stresses involved in the heating process is a further 


possibility. 


There are many difficulties involved in modelling a 
continuous heating process to obtain temperature 
information. First it would be necessary to develop a 
program to solve the thermal conduction problem on the same 
lattice of points as the scatterer, with appropriate 
temperature dependent thermal properties at each point. 
Using the results of this program, a new dielectric data 
lattice would have to be calculated, followed by a new 
lattice of thermal properties, before the process could be 
repeated for the next temperature step. As explained before, 
computer time and storage space are the main limitations to 
this technique. Other mathematical models could no doubt be 
developed, but the resulting set of differential equations 


with variable coefficients would still be difficult to 


solve. 


Discussions with the originator (60) of the modelling 
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technique used in this work suggested that extension to oven 
heating is feasible. Introducing a conductive boundary at 
the edges of the lattice containing the scatterer, or 
providing appropriate lattice truncation conditions would 
Simulate a scatterer within a conducting cavity. Different 
types of wave launching structures could also be simulated 
by assigning appropriate dielectric properties to the points 
of the scatterer in the vicinity of the coupling area. 
Again, computer storage is the most serious limitation to 
application of this model to heating within a microwave oven 


of reasonable size. 


One possible approach towards estimating the stresses 
involved in the heating process is to use the techniques of 
Lin (61) who has compared some of the possible mechanisms of 
microwave auditory effects, as observed by Frey (62). 
According to Lin's calculations, thermal effects and 
strictive forces are the primary causes of the auditory 
effect, and it is reasonable to suppose that these forces 
will be present during microwave heating, where higher 
powers, and therefore higher field strengths, are employed. 
Using the heating model and the appropriate thermal 
properties, it should be possible to calculate the forces 
exerted in different parts of the body. It may then be 
feasible to relate these forces to a threshold for tissue 
damage, which would in turn set an upper limit to the field 


strength permissible in a given application. 
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The results of this work indicate one possible approach 
to the problem of controlled microwave heating. A system 
capable of maintaining a constant heating rate within 
prescribed limits is described and some of the results 
obtained with the system are presented. Various limitations 


of the system and some possible improvements are discussed. 


Measurements on the dielectric properties of solutions 
of water and DMSO are described, and the results obtained 
are presented. A minimum relaxation frequency of 3 Ghz is 
found for a 70% DMSO solution, the relaxation frequency 
increasing to 9 Ghz for 100% DMSO and to 18 Ghz for water. 
The 70% concentration corresponds to the eutectic point of 
the mixture. The minimum relaxation frequency encountered at 
this concentration is apparently due to the oocurrence of 
maximum hydrogen bonding between the water and DMSO 
molecules. In addition, the presence of DMSO at temperatures 
below the freezing point is found to increase the losses of 
the solution, and to reduce the rate of change of €" with 
temperature. The cryobiological implications of these 


results are discussed. 


A numerical model to simulate microwave heating of the 
canine kidney is developed and the model is employed to 


examine the effect of various parameters on heating rate and 
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uniformity. It is found that the addition of DMSO to the 
perfusate enhances the heating rate. The heating uniformity 
is aiso sensitive to the DMSO concentration, and it appears 
that control of the DMSO concentration during heating could 
provide some measure of control of heating uniformity. In 
the vicinity of the renal sinus, which absorbs little 
energy, considerable electric field gradients are found. 
This is due to the low perfusion assumed for this region, 
and suggests that ways should be found to improve perfusion 
in the sinus. The composition of the fluids residing in the 
renal pelvis is aiso found to have considerable effect on 
the heating uniformity. Maximum heating uniformity is found 
to correspond to maximum uniformity of perfusion. Another 
phenomenon which the model revealed is that of strong 
heating of the tubules leading off from the renal pelvis. 
Although the vascular structure is not represented on the 


model, it is supposed that similar effects will occur there. 
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AS explained in Chapter 5, two computer programs were 
used to analyse the raw data obtained from the dielectric 
measurements. A brief description of the function of each 
program is given below, followed by FORTRAN listings of the 


programs. 


alibration Program (CALCOEFF) 


This program takes the values of €" and €" for the 5 
standard solutions, and, using measured values of frequency 
shift and attenuation, performs a bicubic spline fit of the 
data. The spline coefficients are stored for use in the 
second program. Three separate spline fits are performed, 
representing the mean and the upper and lower bounds of the 


standard deviation. 


Using the values of the coefficients obtained in the 
CALCOEFF program, this program calculates €* and €" from the 
frequency shift and attenuation data obtained from 
measurements on the unknown solutions. Again, three values 
are calculated in each case, these being the mean and 


standard deviation bounds. 
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PROGRAM CALCOEFF 


(lini te! el hel (al tal teh a 


QA2a0 


4567 


3456 


5678 


2121 


2323 
2324 


LOCAL OS @  Qi@ 


THIS PROGRAM CALCULATES THE VALUES OF E* AND E" FOR 
UNKNOWN SUBSTANCES. FOR EACH TEMPERATURE, THE MEAN 
AND STANDARD DEVIATION OF NDATA VALUES OF FREQUENCY 
SHIFT AND ATTENUATION ARE FOUND, GIVING DATAL, DATAM 
AND DATAU AS THE ERROR BOUNDS. IN EACH CASE THE 
SUBSCRIPTS ARE USED AS FOLLOWS; 

esse. (I, 1) =TEMPERATURE 

oeee- (1,2) =FREQUENCY SHIFT 

eeee- (1,3) =ATTENUATION 


DIMENSION DATA (100,3) ,DATAL (3) ,DATAM (3) , DATAU (3) 
DIMENSION COEFFL (4) ,COEFFM(4) , COEFFU (4) 
DIMENSION COEPPL (4) ,COEPPM(4) ,COEPPU (4) 
DIMENSION EPL(100),EPM(100) ,EPU(100) 

DIMENSION EPPL (100) , EPPM (100) , EPPU(100) 
LOGICAL*1 ITYPE(4), STRING (16) , TITLE (50) , EQUC, 
BUFFER (60) 


THE COEFFICIENT DATA, AS CALCULATED BY "CALCOEFF", IS 
READ FROM "COEFFICIENTS". 

DO 4567 I=1,5 

CALL MOVEC(12,! ', BUFFER ( ( (I-1) *12) +1) ) 
READ (5,3456) ITLEN, (TITLE(I) ,1=1,50) 
FORMAT (12, 1X,50A1) 

ISTART=(60-ITLEN) /2 

CALL MOVEC (ITLEN, TITLE(1) , BUFFER (ISTART) ) 
WRITE (6,5678) (BUFFER (I) ,1=1,60) 

FORMAT (*1*,60A1) 

READ (5,2121) ITYPE 

FORMAT (4A1) 

IF (EQUC(ITYPE,*TEMP')) GOTO 2323 

CALL MOVEC(16," % CONCENTRATION, STRING) 
GOTO 2324 

CALL MOVEC(16," TEMPERATURE ‘',STRING) 
WRITE (6,1000) (STRING (I) ,I=1, 16) 

READ (8, 1003) (COEFFL (I) -1=1, 4) 

READ (8, 1003) (COEFFM(I) ,1=1,4) 

READ (8, 1003) (COEFFU(I) ,1=1,4) 

READ (8, 1003) (COEPPL (I) ,I=1, 4) 

READ (8, 1003) (COEPPM (I) ,1=1,4) 

READ (8, 1003) (COEPPU (I) ,1=1, 4) 


THE DATA FOR NTEMP DIFFERENT TEMPERATURES (IN ORDER 
OF INCREASING TEMPERATURE) IS NOW READ IN. THERE ARE 
NDATA FREQUENCY SHIFTS AND ATTENUATIONS FOR EACH 


TEMPERATURE. 
MEANS AND STANDARD DEVIATIONS ARE THEN CALCULATED 


FOR EACH TEMPERATURE 


READ (5,1002) NDATA,NTEMP 
DO 130 L=1,NTEMP 
SUMF=0 
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100 


AaAAN 


120 


130 


SUMA=0 

SUM2F=0 

SUM2A=0 

DO 100 I=1,NDATA 

READ (5, 1004) (DATA (L,J) ,J=1, 3) 
SUMF=DATA (L,2) +SUMF 
SUMA=DATA (L, 3) +SUMA 
SUM2F=DATA (L, 2) **2.+SUM2F 
SUM2A=DATA (L, 3) **2.+SUM2A 
CONTINUE 

AK=SUMF*SUMF/NDATA 

BK= (ABS (SUM2F-AK) ) /(NDATA-1) 
SDF=SQRT (BK) 
AA=SUMA*SUMA/NDATA 

BB=(ABS (SUM2A-AA) ) /(NDATA-1) 
SDA=SQRT (BB) 

DATAL (2) =SUMF/NDATA-SDF 
DATAL (3) =SUMA/NDATA-SDA 
DATAM (2) =SUMF/NDATA 

DATAM (3) =SUMA/NDATA 

DATAU (2) =SUMF/NDATA+SDF 
DATAU (3) =SUMA/NDATA+SDA 


VALUES OF E" AND E™ AND THEIR LOWER AND UPPER 
BOUNDS ARE NOW CALCULATED FOR EACH TEMPERATURE. 


K1=K-1 

EPL (L) =0 

EPM(L) =0 

EPU (L) =0 

EPPL (L) =0 

EPPM (L) =0 

EP PU (L) =0 

DO 110 I=1,4 

EPL (L) =EPL (L) +COEFFL (I) *DATAL (2) ** (I-1) 

EPM(L) =EPM(L) +COEFFM (I) *DATAM (2) **(I-1) 

EPU (L) =EPU (L) +COEFFU (I) *DATAU (2) **(I-1) 

DO 120 I=1,4 

EPPL (L) =EPPL (L) +COEPPL (I) *ALOG 10 (EPL (L-1) ** (I-1) 
EPPM (L) =EPPM (L) *COEPPM (I) *ALOG 10 (EPM (L-1) ** (I-1) 
EPPU (L) =EPPU (L) +COEPPU (I) *ALOG 10 (EPU (L-1) ** (I-1) 
EPPL (L) =10. ** (DATAL (3) /20+EPPL (L-1. 25) 

EPPM (L) =10. ** (DATAM (3) /20+EPPM (L-1. 25) 

EPPU (L) =10. ** (DATAU (3) /20+EPPU (L-1. 25) 
TANL=EPPL (L) /EPU (L) 

TANM=EPPM (L) /EPM(L) 

TANU=EPPU (L) /EPL (L) 

PRINT 1005, EPL (L) ,EPPL(L) , TANL 

PRINT 1006, DATA(L, 1) ,EPM(L) ,EPPM(L) , TANM 

PRINT 1007, EPU (L), EPPU(L) ,TANU 
WRITE DATA FOR PLOT ROUTINE 

WRITE (3, 2000) DATA(L,1),EPL(L) ,EPM(L), EPU (L) 
WRITE (4,2000) DATA(L,1), EPPL(L), EPPM(L), EPPU(L) 
CONTINUE 
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FORMAT ("0',T8, 16A1,1727,"EP*,T41,"EPP",1T53,"LOSS TAN") 
FORMAT(* *,2F7.2) 

FORMAT (11,14) 

FORMAT (1X,6E15.4) 

FORMAT (3F8. 1) 
FORMAT(*0',T25,F5.2,T40,F5. 2,754, F6. 3) 

FORMAT(* *,1712,F6.1,1725,F5.2,T40,F5.2,T54,F6.3) 
FORMAT(* *,125,F5.2,T40,F5.2,T54, F6.3) 

FORMAT(* ',3F8.1) 

FORMAT (4F6. 2) 

STOP 

END 
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PROGRAM CALCULATE 


THIS PROGRAM CALCULATES THE COEFFICIENTS FOR A LEAST 
SQUARES FIT OF FREQUENCY SHIFT DATA (FS) TO E* (EP). 
3 DISTINCT CURVES (MEAN, MEAN+S.DEVIATION AND 
MEAN-S.DEVIATION FOR N SEPARATE READINGS) ARE FOUND, 
ENABLING THE STANDARD DEVIATION OF ANY 

MEASUREMENT ON UNKNOWN SUBSTANCES TO BE ESTIMATED. 


DIMENSION EP(7),EPP(7),FS(7,10),AT(7, 10) 
DIMENSION FU(7) ,FL(7) ,FM(7) 

DIMENSION AL(7) ,AM(7) ,AU(7) 

DIMENSION COEFFL (4) ,COEFFM(4) , COEFFU (4) 
DIMENSION C(6),S(6),A(3) ,B(3) ,X(7) 

DIMENSION EP1(7),EPP25L(7) , EPP25M (7) , EPP25U (7) 
DOUBLE PRECISION P(14),T (16) 


N (MAX. 20) IS THE NUMBER OF READINGS FOR EACH 
CALIBRATION LIQUID AND NCL (MAX. 10) IS THE NUMBER 
OF CALIBRATION LIQUIDS. 


N=6 

NCL=7 

READ 1001, (EP (Z) ,1=1, NCL) 
READ 1001, (EPP (I) ,I=1,NCL) 
DO 100 I=1,NCL 

READ 1002, (FS (I,J) ,J=1,N) 
READ 1002, (AT (I,J) ,J=1,N) 
CONTINUE 

DO 120 I=1,NCL 

FS1=0 

FS2=0 

AT1=0 

AT2=0 


CALCULATE THE MEANS AND STANDARD DEVIATIONS OF THE 
INPUT DATA AND PLACE THEM IN THE APPROPRIATE 
COEFFICIENT MATRIX. 


DO 110 J=1,N 
FS1=FS1+FS (I,J) 
FS2=FS2+FS (I,J) **2. 
AT1=AT1+tAT (I,J) 
AT2=AT2+AT (I,J) **2. 
CONTINUE 
AK=FS1*FS1/N 
BK=(FS2-AK) /(N-1) 
SDF=SQRT (BK) 
AA=AT1*AT1/N 

BB= (AT2-AA) / (N-1) 
SDA=SQRT (BB) 

FU (1) =FS1/N+SDF 

FM (I) =FS1/N 

FL (I) =FS1/N-SDF 


} 
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AU (I) =AT1/N+SDA 
AM (I) =AT1/N 

AL (I) =AT1/N-SDA 

CONTINUE 

DO 260 J=1,3 

IF (J-2) 130,150,170 

DO 140 I=1,NCL 

X (I) =FU (I) 

GO TO 190 

DO 160 I=1,NCL 

X (I) =FM (I) 

GO TO 190 

DO 180 I=1,NCL 

X (I) =FL (1) 

CONTINUE 

CALL RLFOTH (X,EP,NCL,100.,3,1D,P,C,S,A,B,IER) 
CALL RLDOPM (C,ID,A,B,T) 
KD=ID+1 

IF (J-2) 200,220,240 

WRITE (6, 1004) (C(I) ,1=1,KD) 
GO TO 260 

WRITE (6, 1004) (C(I) ,2=1,KD) 
GO TO 260 

WRITE (6, 1004) (C(I) ,I=1,KD) 
CONTINUE 


ESTABLISH VECTORS OF LOG(E"-1) (EP1) AND LOG(E") 
AT 25 DB (EPP25). 


DO 270 I=1, NCL 

EP 1(1I) =ALOG10 (EP (I-1) 

EPP25L (I) =ALOG 10 (EPP (I) -AU (I) /20+1. 25 
EPP25M (I) =ALOG10 (EPP (I) AM (I) /20+1. 25 
EPP25U (I) =ALOG 10 (EPP (I) -AL (I) /20+1. 25 
CONTINUE 


CALCULATE THE LEAST SQUARES COEFFICIENTS FOR THE ABOVE 
VECTORS FOR THE THREE CASES, MEAN-SD, MEAN AND MEAN+tSD 
COEPPL,COEPPM AND COEPPU ARE THE COEFFICIENT VECTORS 
FOR THE CURVES REPRESENTING THE MEAN-SD, MEAN AND 
MEAN+SD OF EP-1 VS. EPP25 


DO 280 I=1,NCL 
X (I) =EP1(I) 

CALL RLFOTH (X,EPP25L,NCL, 100.,3,1D,P,C,S,A,B,IER) 
CALL RLDOPM (C,ID,A,B,T) 

KD=ID+1 

PRINT 1004, (C(I) ,1=1,KD) 

DO 290 I=1,NCL 

X (I) =EP1 (I) 

CALL RLFOTH (X,EPP25M,NCL,100.,3,1D,P,C,S,A,B,1IER) 
CALL RLDOPM (C,ID,A,B,T) 

KD=ID+1 

PRINT 1004, (C(I) ,1=1, KD) 
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DO 300 I=1,NCL 

X (I) =EP1 (I) 

CALL RLFOTH (X,EPP25U, NCL, 100.43, 1D, PeC,yS,AyB,1ER) 

CALL RLDOPM (C,ID,A,B,T) 

KD=ID+1 

PRINT 1004, (C(I) ,I=1, KD) 

FORMAT (11) 

FORMAT (F5.2,1X,F5.2, 1Xg F502, 1X,F 5-29 1X, F522, 1X, 
CF5.2, 1X, F5.2) 

FORMAT (F5.2,1X,F5.2, 1X9 F 5029 1X,F502, 1X, F522, 1X, F5~ 2) 

FORMAT ("O",*IER =",13) 

FORMAT(* ',6515.4) 

FORMAT(*O", "VALUES OF EP*,7F6. 2) 

FORMAT(*O*," VALUES OF EPP*,7F6.2) 

FORMAT(*O*, "VALUES OF FS‘) 

FORMAT(*0*,6F6.1) 

FORMAT(*0",* VALUES OF AT*) 

FORMAT (*0%, *FU (",11,")=",F5.1,5X,"FM("g11,") =", F501, 
C5X,"FL(*, I1,°)=",F5.1) 

FORMAT(*0*, *AU(*, 11,*) =",F5.1,5X, "AM(",I1,")=",F5.1, 
CSP! NEU 117" yHt7F5e 1) 

STOP 

END 
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Appendix B 


As explained in Chapter 8, a program to calculate the 
heating index of a plane incident microwave field on an 
arbitrary shaped body was written based on a similar program 
due to faflove (T3). This program is called SCAT.KID.2 and a 


FORTRAN listing is given below. 


Input to the program includes the frequency, the total 
humber of steps for which the program is to be run, the 
humber of steps for which the envelope is to be scanned and 
the step ae which scanning is to begin. Also included are 
the number of media and their dielectric properties, and a 
media lattice describing the distribution of the various 
media. Output from the program is an index fron 0-9 
crepresenting the amount of energy dissipated at each point 
within the lattice. The maximum dissipated power in W/m, 


corresponding to an index of 9, is also returned. 
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PROGRAM SCAT.KID.2 


——— eS Se Se SS Se 


QQ’ Qa) O.@@ QQ QaaGiai@e 


1110 
1103 


Cc 
C 


1101 


1104 


DAWDOLEWH | 


PROGRAM SCAT.KID.2 


STEADY TEM IRRADIATION OF A CANINE KIDNEY 

INCIDENT WAVE HAS THE COMPONENTS EZ,HX,KY 

5 DIFFERENT DIELECTRICS ARE DEFINED 

27 X 34 X 15 CELL YEE LATTICE IS USED 

UNIT CELL DIAMETER = 0.05*LAMBDA (MIN) =DX 

EVEN SYMMETRY ABOUT PLANES X=27.5*DX AND Z=15.0*DX 
IS ASSUMED AND ADAPTIVE SOFT LATTICE TRUNCATIONS 
ARE USED. 


THE MEDIA LATTICE AND THE ENVELOPE PRINTOUT 
ARE TIED TO UNIT 6. 
CONTROL INPUT IS ON UNIT 5, DATA INPUT ON 8. 


REAL MUZ 

COMPLEX JKAY, CMPLX, CSORT 

INTEGER IPR(27),N(27) ,IVF(9) 

REAL EX0(15,27), EZ0(15,27), EXT1F(15,27), 
EXT2F (15,27), EYT(34,14,27) , ENY (34, 15,27), 
EZT1F (15,27), EZT2F (15,27), EXT(34, 14,27), 
HYT (34,14,15), HZT(34, 14,15), ENX(34,15,27), 
ENZ(34,15,27), EX(34,16,28), EY(34,16,28), 
HX (34,16,28), HY (34,16,28), HZ(34, 16,28), 
EZ (34,16,28), MEDIA(34,15,27), EPS(9), 
EXTR(14,15,27), EZTR(14, 15, 27) 
SIG (9) ,CA (9) ,CB(9) 


oes--I. PROBLEM PARAMETERS..... 
READ (5,1100) FREQ 

FORMAT (E6.1) 

PRINT 1102, FREQ 

FORMAT(*1*,"*THE FREQUENCY IS = ',E9.2,* HZ") 
DX=1.5E-3 


MAXDT IS THE TOTAL NO OF STEPS REQUIRED. 
MINDT IS THE TIME AT WHICH ENVELOPE SCANNING BEGINS. 
MIDDT IS THE NUMBER OF STEPS FOR WHICH THE ENVELOPE 
IS SCANNED. 

READ (5,1110) MAXDT, MINDT, MIDDT 
FORMAT (13, 1X,13, 1X,13) 

PRINT 1103, MAXDT, MINDT, MIDDT 

FORMAT(* *,*MAXDT = ',13,/,"MINDT = ',13, 
1/*MIDDT = ',12) 


S85, SAT Te EDTELECTRIC MEDIA... 
MPR=5 

READ (5,1101) (EPS (I) ,1=1,MPR) 
READ (5,1101) (SIG(I) ,1=1,MPR) 
FORMAT (9 (F5.2, 1X) ) 

DO 1104 I=1,MPR 

PRINT 1105, I, EPS (I) /1,SiG (I) 
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1105 FORMAT (! Dre nha (ML te y= itr 5 2. 5x 


EQUATE MTS AND FORTRAN LOGICAL UNITS 
CALL FINCMD(*EQUATE 9=9*,10) 


++e--III. BASIC AND DERIVED CONSTANTS..... 
PI=3.1415927 
MUZ=4.0*PI*1.0E-7 
EPSZ=8. 854E-12 
OMEGA=2.0*PI*F REQ 
DT=DX/5.99585E+8 
MDT=MINDT+MIDDT 
R=DT/EPSZ 
RA=DT*DT/DX/DX/MUZ/EPSZ 
RB=DT/DX/MUZ 
RC=MUZ*EPSZ*OMEGA**2 
RD=MUZ*OMEGA 
RE=5.99585E+8/0MEGA 
RF=1000.0*RB 
RG=OMEGA*DT 
RH=2. O#RF**2 
FMPR=FLOAT (MPR) 
NEEST=0 
DO 1 IAA=1,27 

1 N(IAA)=IAA 


eeeeeIV. CALCULATE THE PROPAGATION CONSTANT AND..... 
HENCE THE MAXIMUM DELAY FOR EACH TISSUE TYPE 
DO 2 IB=1,MPR 
JKAY=CMPLX (0.0, 1.0) *CSQRT (CMPLX (EPS (IB) *RC, 
1 SIG (IB) *RD) ) 
IVF (IB) =AMIN1 (REXAIMAG (JKAY-0.5, 13.5) 
CA (IB) =1. 0-R*SIG (IB) /EPS (IB) 
2 CB(IB)=RA/EPS (1B) 


«ee..V. LOAD THE STORAGE VECTOR..... 
ese.eA. ZERO INITIAL FIELDS (FRONT REAR)...... 
DO WSAICR=1727 
pO 3 Ic=1,15 
EXO (IC, ICA) =0.0 
3 EZO(IC,ICA) =0.0 
DONIO10pLECA=1, 27 
DO 3010 Ic=1,15 
EXT1F (IC, ICA) =0.0 
3010 EZTIF (IC,ICA) =0.0 
DO 3011 ICA=1, 27 
DO S01181C21/15 
EXT2F (IC, ICA) =0.0 
3011 EZT2F (IC,ICA) =0.0 
DO 301: ICB=1,27 
DO 301 IcA=1,15 
pO 301 Ic=1,14 
EXTR (IC, 1CA,ICB) =0.0 
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EZTR (IC, ICA, ICB) =0.0 


«+---B. TYPE OF DIELECTRIC MEDIUM..... 
DO 4 IK=1,15 

DO 4 ID=1,34 

IE=16-1K 

READ (8, 1113) (MEDIA (ID, IE, IEA) , LEA=1, 27) 
FORMAT (27F1.0) 

CONTINUE 


eoeeeC. ZERO INITIAL FIELDS (CENTRE)..... 
DO 402 IC=1,27 

DO 402 IB=1,14 

DO 402 IaA=1,34 
EXT(IA,1IB,IC) =0.0 
EYT(IA,1IB,IC) =0.0 
DO 4020 Ic=1,15 
DO 4020 IB=1,14 
DO 4020 IA=1,34 
HYT (IA, 1B,1C) =0.0 
HZT(IA,1IB,IC) =0.0 
DO 403 Ic=1,27 

DO 403 IB=1,15 

DO 403 IA=1,34 
ENX(IA,1IB,IC) =0.0 
ENY (IA, IB,IC) =0.0 
ENZ(IA,1B,IC) =0.0 
DO 4030 IC=1,28 
DO 4030 IB=1,16 
DO 4030 IA=1,34 
EX (IA, IB, IC) =0.0 
EY (IA, IB, IC) =0.0 
EZ (ZA, IR, 1c) =0-0 
DO 4031 IC=1,28 
DO 4031 IB=1,16 
DO 4031 IA=1,34 
HX (IA, 1B, IC) =0.0 
HY (IA, 1B, IC) =0.0 
HZ (IA, 1B, IC) =0.0 


seeeeVIe PRINT THE DIELECTRIC MEDIA LATTICE..... 

DO 8 IF=1,15 

PRINT 6,1F 

FORMAT (1H1,////////455X-"DIELECTRIC MEDIA LATTICE', 
J fPGQLIYBUANEOKE= 1,1 27//7,38X, 10" 7/7) 

DO 67 IG=1,34 

IH=35-IG 

DO 66 II=1,27 

IPR(Z1I)= MEDIA (1IH,IF,11I) 

IF (IPR(II).LT.1.0.OR.IPR(II).GT. FMPR)GO TO 70 

CONTINUE 

PRINT 7,1H, (IPR (1d) ,1J=1, 27) 

FORMAT (37X,12,4X,2712) 

PRINT 9, (N(IK) pLK=1,9)» (N(IL) - IL=1, 9)» (N (IM) , IM=1, 7) 
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9 FORMAT (///e@4O0X, "I", 2X,912,10(" 1%) ,8(8 2) G7 
161X," 0*,912," 0*,722) 


eeee VIL. TIME STEPPING LOOP..... 
DO 65 J=1,MAXDT 


eee-e-e ESTABLISH LOOP AND FIELD PARAMETERS..... 


109° TERM=RF*SIN (FLOAT (J) *RG) 


MAXY=MINO (6+IFIX (FLOAT (J) /2.0) ,34) 


ecee eA. E FIELD LOOP..... 
DO 22 JA=1,MAXY 
IF (JA.EQ.34)GO TO 14 


wae d. Vo EXSETERATIONGS 29. 
IF (J.GE.MINDT) GO TO 111 
DO 110 JCA=2,27 
DO 110 JBA=2,15 
NF=MEDIA (JA, JBA, JCA) 
110 EX (JA,JBA, JCA) =CA (NF) *EX (JAJJBA, JCA) +CB (NF) * 
1 (HY (JA,JBA-1,JCA) - HY(JA,JBA, JCA) + 
2 HZ (JA+1,JBA,JCA-HZ(JA, JBA,JCA)) 
GO TO 112 
111 DO 10 JC=2,27 
DO 10 JB=2,15 
NF=MEDIA (JA,JB,JC) 
EX (JA, JB, JC) =CA (NF) *EX (JA, JB, JC) +CB (NF) * 
1 (HY (JA,JB-1,3C) - HY(JA,JB,JC) + 
2 HZ (JA+1,JB, JC-HZ (JA, JB,JC)) 
10 ENX (JA, JB, JC) =AMAX1 (ABS (EX (JA, JB, JC) ) ,ENX (JA, JB, JC) 
112 REX=EX (JA,2,10) /RB 
IF (ABS (REX) -GT.1.0E+4)GO TO 42 


eeeee2. EX SOFT LATTICE TRUNCATION..... 
EX (JA, 1,2) = (EXT (JA, 1, 2) #EXT (JA, 1,3)) /2.0 
EX (JA, 1,27) =(EXT (JA, 1,26) +EXT (JA, 1,27)) /2.0 
DO 11 JD=3,26 
11 EX (JA, 1, 5D) = (EXT (JA, 1,3D-1) +EXT(JA, 1,JD) 
1 +EXT(JA,1,5D+1) ) /3.0 
DO 13 JE=2,27 
NF=MEDIA (JA, 1, JE) 
JF=IVF (NF) 
DORI2eaG=1,0F 
12 EXT(JA,JG, JE) =EXT (JA, JG+1,JE) 
13 EXT(JA, JG+1, JE) =EX (JA,2,JE) 


aes 235 VEG ITERATION. (00 
14 IF(J.GE.MINDT) GO TO 116 
DO 115 JIA=1,27 
DO 115 JHA=2,15 | 
NF=MEDIA (JA,JHA,JIA 
BY (92, JHA, FEA) =CA (NF) *EY (JA, HA, JIA) +CB (NF) * 

(HX (JA, JHA, JIA-HX (JA, JHA-1, JIA) 

2 +HZ (JA, JHA, JIA-HZ (JA, JHA, JIA+1) ) 
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115 CONTINUE 
GO TO 117 
116 DO 15 JI=1,27 
DO 15 JH=2,15 
NF=MEDIA (JA,JH, JI) 
EY (JA, JH, JL) =CA (NF) *EY (JA, JH, JI) +CB (NF) * 
1 (HX (JA, JH, JI-~HX (JA, JH-1, JI) 
2 +HZ(JA,JH, JI-HZ(JA, JH, JI+1)) 
15 ENY(JA,JH, JI) =AMAX1 (ABS (EY (JA, JH, JI)) , ENY (JA, JH, JI) ) 
117 REY=EY (JA,2,10) /RB 
IF (ABS (REY) .GT.1.0E#4)GO TO 42 


oeeee4. EY SOFT LATTICE TRUNCATION..... 
EY (JA, 1, 1) = (EYT (JA, 1,1) +EYT (JA, 1, 2) ) /2.0 
EY (JA, 1,27) =(EYT (JA, 1,26) +2. O*EYT (JA, 1, 27)) /320 
DO 16 JJ=2,26 
16 EY (JA,1,Jd)=(EYT(JA,1,JJ-1) +EYT(JA, 1,Jd) 
1 +EYT(JA,1,J3+1)) /3.0 
DO 18 JK=1,27 
NF=MEDIA (JA, 1, JK) 
JF=IVF (NF) 
DO 17. JL=1,JF 
17 EYT(JA, JL, JK) =EYT (JA, JL+1, JK) 
18 EYT(JA,JL+1,dK) =EY (JA,2,JK) 
IF (JA.EQ.34)GO TO 23 


£98595 EL METERATIONL SS .2 
IF (J.GE.MINDT) GO TO 120 
DO 119 JNA=1,27 
DO 119 JMA=1,15 
NF=MEDIA (JA, JMA, JNA) 
EZ (JA, JMA, JNA) =CA (NF) *EZ (JA, JMA, JNA) +CB (NF) 
1 * (HX (JA, JMA, JNA-HX (JA+1,JMA,JNA) 
2 +HY (JA, JMA,JNA+1-HY (JA, JMA,JNA)) 
119 CONTINUE 
GO TO 121 
120 DO 19 JN=1,27 
DO 19 JM=1,15 
NF=MEDIA (JA,JM,JN) 
EZ (JA, JM, JN) =CA (NF) *EZ (JA, JM, JN) +CB (NF) * (HX (JA, JM, JN) 
1 -HX (JA+1,3M,JN) tHY (JA, JM, JN+1) 
2 -HY (JA,JM,JN) ) 
19 ENZ(JA,JM, JN) =AMAX1 (ABS (EZ (JA, JM, JN) ) ENZ (JA, JM, JN) ) 
121 REZ=EZ(JA, 1,10) /RB 
IF (ABS (REZ) .GT.1.0E+4)GO TO 42 
IF (JA-EQ.1)GO TO 122 
LIF (JA-3) 22, 123,22 


.-cee6. EX,EZ TRUNCATIONS AT FRONT LATTICE FACE..... 
122 DO 1221 JzZB=1,27 
DOMIoo eo ZA-1710 
EXO (JZA,JZB) =EXT1F (JZA,J5ZB) 
EXT1F (JZA,JZB) =EXT2F (JZA,JZB) 
1221 EXT2F (JZA,JZB) =EX (1,5ZA,JZB) 
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DO 1222 JZB=1,27 

DO 1222 Jza=1,15 
EZ0(JSZA,JZB) =EZT1F (JZA,JZB) 
EZT1F (JZA,JZB) =EZT2F (JZA, JZB) 
EZT2F (JZA,JZB) =EZ (1,3 ZA, J3ZB) 
GO TO 21 


eese+7. SOFT TEM WAVE SOURCE CONDITION..... 

DO 20 JP=1,27 

DO 20 JO=1,15 

EZ (JA, JO, JP) =TERM+EZ (JA,JO, JP) 

CONTINUE 

CONTINUE 

GO TO 27 

CONTINUE 


eeeee8. EX, EZ TRUNCATIONS AT REAR LATTICE FACE..... 


DO 230 JSA=1,27 

DO 230 JRA=1,15 

EX (34, JRA, JSA) =EXTR(1,JRA,JSA) 
EZ (34, JRA, JSA) =EZTR(1,JRA,JSA) 
DO 25. JS=1,27 

DO 25 JR=1,15 
NF=MEDIA (34, JR,JS) 

JF=IVF (NF) 

DO 24 JT=1,JF 

EXTR (JT, JR, JS) =EXTR(JT+1,3R,JS) 
EZTR (JT, JR, JS) =EZTR(JT+1,5R,JS) 
EXTR (JF#1,JR, JS) =EX (33,3R,JS) 
EZTR (JF+1,3R, JS) =EZ(33,J3R,dS) 


joes DSRATFIELDILOOPT 4.2 
DO 4O K=1, MAXY 
IF(K.GE.2)GO TO 28 
GO TO 29 


f.Sbet. \HX ITERATION... ..< 
DO 30 KB=1,27 
DO 30 KA=1,15 
HX (K, KA, KB) =HX (K, KA, KB) +EY (K, KA+ 1,KB-EY (K, KA, KB) 
+EZ (K-1,KA,KB-EZ (K, KA, KB) 


30 CONTINUE 


GOTO 306 


29 DO 305 KBZ=1,27 


1 


DO 305 KAz=1,15 
HX (K, KAZ, KBZ) =HX (K, KAZ, KBZ) +EY (K,KAZ+1, KBZ) 
- =BY (K, KAZ, KBZ) +EZO (KAZ, KBZ-EZ (K, KAZ, KBZ) 


305 CONTINUE 
306 CONTINUE 


IF (K.EQ.34)GO TO 35 


eeeee2- HY ITERATION...~-«- 
DO 31 KD=2,27 
DO 31 KC=1,15 
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HY (K, KC, KD) =HY (K, KC, KD) +EX (Ky KC, KD~ EX (K, KC+#1, KD) 
1 +EZ (K, KC, KD-EZ (K, KC, KD-1) 
31 CONTINUE 


ooe+-3. HY SOFT LATTICE TRUNCATION..... 
HY (K,1,1)=(HYT (K, 1,1) tHYT(K, 1,2) ) /2.0 
HY (K, 15,1) =(HYT(K, 1,14) +2. 0*HYT(K,1, 15) ) /3-0 
DO 32 KE=2,14 
32 HY (K, KE, 1)=(HYT(K, 1, KE-1) +HYT (K, 1, KE) 
1 +HYT (K,1,KE+1)) /3.0 
DO 34 KF=1,15 
NF=MEDIA (K,KF, 1) 
KG=IVF (NF) 
DO 33 KH=1,KG 
33 HYT(K,KH,KF) =HYT(K,KH+1, KF) 
34 HYT(K, KH+1,KF) =HY (K,KF,2) 


va oe AUNGHZ OFTERATIONSOUS. 

IF (K.GE.2)GO TO 35 
DO 341 KJ=2,27 
DO 341 KI=2,15 

341 HZ (K,KI, KJ) =HZ (K, KI, KJ) + EX (Ky KI, KJ- EXO (KI, KJ) 
1 +EY (K, KI, KJ-1-EY (K, KI, KJ) 
GO TO 351 

35 DO 36 KJ=2,27 
DO 36 KI=2,15 

36 HZ (K,KI,KJ) =HZ (K, KI, KJ) +EX (K, KI, KJ-EX (K-1, KI, KJ) 
1 +EY (K, KI, KJ-1-EY (K, KI, Kd) 

351 CONTINUE 


wees-5- HZ SOFT LATTICE TRUNCATION..... 
HZ (Ky 2,9 1) =(HZT (K, 1,2) +HZT (K, 1,3) ) /2.0 
HZ (K,15, 1) = (HZT (K, 1,14) +HZT (K, 1,15) ) /2.0 
DO 37 KK=3,14 
37 HZ (K,KK, 1) = (HZT (K, 1, KK-1) #HZT (K, 1, KK) 
1 +HZT(K, 1, KK+1)) /3-0 
DO 39 KL=2,15 
NF=MEDIA (K, KL, 1) 
KG=I1 VF (NF) 
DO 38 KM=1,KG 
38 HZT(K, KM, KL) =HZT (K, KM+1,KL) 
39 HZT(K, KM+1,KL) =HZ (K, KL, 2) 
40 CONTINUE 


eeeeCe PRINT-OUT ROUTINE..... 
DO 41 L=MDT,MAXDT,MIDDT 
IF (J.EQ.L)GO TO 43 
41 CONTINUE 
GO TO 64 
42 NTEST=1 
FIND THE NORMALISING PARAMETER..... 


ra Re 
43 EMAX=0.0 
pO 45 LC=1,27 
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DON5°iE=4,15 
DO 45 LA=1,34 
LD=MEDIA (LA, LB, LC) 
ER= (ENX (LA, LB, LC) **2+ENY (LA,LB,LC) **2 
1+ENZ(LA,LB,LC) **2) *SIG (LD) /RH 
IF (ER-EMAX) 45,45, uu 
44 EMAX=ER 
45 CONTINUE 
PRINT 46, EMAX 
46 FORMAT('1*,"THE MAXIMUM DISSIPATED POWER IS ',EF10.3, 
11X,"WATTS/M**3 FOR UNIT INCIDENT ELECTRIC FIELD‘) 
EMAX=EMAX/9 
PRINT 47, EMAX 
47 FORMAT(* *,*THE NORMALISING FACTOR IS ',E10.3) 


eoe--2. CALCULATE THE DISSIPATED POWER AT BACH ..... 
POINT, NORMALISE IT TO AN INDEX IN THE 
RANGE 0-9, AND PRINT OUT THE INDEX IN 
THE SAME FORM AS IN THE MEDIA LATTICE. 
DO 52 LF=1,15 
PRINT 48, J, LF 
48 FORMAT (1H1,////////,46X," NORMALISED HEATING INDEX AT', 
1* TIME STEP ',13,//,62X,"PLANE K = ',12,//,38X,'d3',//) 
DO 50 LH=1,34 
LE=35-LH 
DO 49 LG=1,27 
LI=MEDIA (LE,LF, LG) 
49 IPR(LG) = (ENX (LE, LF, LG) **2+ENY (LE, LF, LG) **2 
1+ENZ (LE, LF, LG) **2) *SIG (LI) / (RH*EMAX) 
50 PRINT 51,LE, (IPR(LJ) , LJ=1, 27) 
51 FORMAT (37X,12, 4X, 2712) 
52 PRINT 53, (N(LK) ,LK=1,9),» (N(LL) ,LL=1,9) , (N (LM) ,LM=1,7) 
53 FORMAT (///,40X,*1",2X,912,10(* 1%),8(* 2%),/,61X," Of, 
1912," 0*,712) 


wage 5. DEFAULT OPTION..... 
IF (NTEST.EQ.1)GO TO 70 


seece6e ENVELOPE RESET TO ZERO..... 

DO 62 LV=1,27 

DO 62 LU=1,15 

DO 62 LT=1,34 

ENX (LT, LU, LV) =0.0 

ENY (LT, LU, LV) =0.0 
62 ENZ(LT,LU,LV) =0.0 
64 CONTINUE 
65 CONTINUE 
70 STOP 

END 
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Appendix C 


selected Physical Properties _of DMSO 
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ee 8 
l Property | Value | 
{Boiling Point at 760mm Hg, °C 1189.0 | 
[Vapour pressure at 20°C, mm Hg {0.37,0.417 | 
ISpecific heat at 25°C, cal/gn 1.4698 | 
[Specific heat of vaporization at b.p., cal/gm {|175 | 
[Molar Heat of vaporization at b.p., kcal/mole {13.67 | 
{Entropy of vaporization at b.p., cal/deg.mole |29.6 | 
[Enthalpy of vaporization at 25°C, kcal/mole {12.64 | 
|Coefficient of expansion, ml/gn.deg {0.00088 i 
[Specific gravity at 25°C, g/ml 11.096 | 
[Melting point, °C {18.55 | 
|Specific heat of fusion, cal/gm {38.8 | 
{Molar heat of fusion, kcal/mole 13543 ,3.03 >} 
{Entropy of fusion, cal/deg.mole {10.4,11.7 | 
{Specific heat of solid at 18.45°C, cal/gn }0.5 { 
{Heat of formation, kcal/mole {-47.7 | 
|Heat of combustion, cal/gn 16050 | 
{Flash point, open vessel, °C 195 | 
{Molai £.p. depression, deg/mole/1000gn 14.4,4.36 | 
jRefractive index (nD2@®) {1.4783 | 
jRefractive index (nD25) {1.47674 | 
{Molar refractivity, Rm $20.12 | 
{Polarizability, cm3 {7.97x10-2¢ | 
{Dielectric constant at 20°C (€static) 148.9 | 
[Dielectric constant at 25°C (€static) [46.4 | 
{Dipole moment, 25°C, D {4.11 | 
| Viscosity, 20°C, cP 12.473 | 
{Viscosity, 25°C, cP {1.99 | 
{Surface tension at 20°C, dyne/cm 146.2 | 
jSurface tension at 25°C, dyne/cn {42.86 | 
{Specific conductance at 20°C, ohm-1!cm—? {3x10-8 | 
E npr nea S 


The above properties are obtained from (S4), page 5. 


This source lists the references in which the above data was 


pubiished. 
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Appendix D 
In order to make available the numerical data obtained 
from the dielectric measurements, the tables in this section 
are presented. 


some comments regarding their presentation 


are in order. 


Fourteen different sets of measurement parameters were 


used, as outlined below: 
(1) Distilled water 
(2) 1% DMSO 

{3) 5% DMSO measurements in the 
(4) 10% DMSO temperature range 
(5) 20% DMSO -60°C to +100°C 
(6) 70% DMSO 


(7) 100% DMSO 


(10) 1.0% Nacl + 10% DMSO -60°C to +40°C 
(11) 5% dextran + 10% DMSO ea, = Sa ee a 
“eoxpravee hee concentrations in 10% 
(13) DMSO at 23°C steps from 0% to 100% 
(14) DMSO at 50°C at each temperature 
The data obtained in the measurements has been reduced 


5% dextran 


0.1% NaCl + 10% DMSO 
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measurements in the 


temperature range 
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to values of €", €", loss tangent and conductivity (mho). 
The standard deviation as a percentage of the value is 
quoted in the case of €' and €", except for cases where this 
exceeds 100% as in the vicinity of the freezing point of 
some mixtures. In these cases, the deviation is left blank. 
Standard deviations are not quoted for the derived 
quantities, loss tangent and conductivity, because they are 
not measured directly. The quoted value is calculated from 
the mean values of €" and €". In the case of sets 1-11, the 
figures quoted for decreasing temperatures were those 
obtained during cooling, and the figures quoted for 
increasing temperatures were those obtained during heating. 
The differences are attributable to supercooling effects, 
and to the existence of mixtures of solid and liquid 
material just beyond the state transition point. For 
purposes of plotting the curves of chapter 6, some liberties 


were taken in these regions, as described in that chapter. 
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DIELECTRIC DATA FOR DISTILLED WATER 


TEMP. 
DEG.C 


24.5 
21.5 
Woes 
10.5 
15 
ci U 
me OILY 
= ES 
ae 
-15.0 
Ueno 
=43.5 


EP 


76.58+t- 1.9% 
79.04+- 1.8% 
80.20+- 2.0% 
81.09+—- 1.9% 
81.23+- 1.8% 
77.90+- 4.0% 
65.72+- 7.2% 
41.85+-13. 3% 
18. 56+-18. 2% 
4.21+-29.7% 
3.114+-10. 9% 
3.04+4+-12. 5% 
3.06+-12.4% 
2.924-12. 7% 
2.664+-14.3% 
2.80+-13. 2% 
2.64+-13. 6% 
2- 804-14. 3% 
2.-834-14. 5% 
2.80+-16. 1% 
2.80+-16. 1% 
2.87+-13. 6% 
4.394-31.4% 
8.99+-26.8% 
15.46+-18.5% 
22. 76+-16.5% 
30. 86+-15.9% 
41.014-14. 3% 
52.03+-11.4% 
63.12+- 9.1% 
720-25t- 5.3% 
76.45+- 5.6% 
78.U44+- 4.7% 
77.97+- 4.7% 
75.22+- 2.5% 
73.69t- 2.8% 
72.31+- 4.5% 
71.38+- 4.2% 
70.65t+- 4.2% 
68.97+- 4.5% 
65.50+- 2.5% 
64.43+- 2.3% 
63-39+- 2.6% 
62-09t- 2.7% 
61.40+- 2.5% 
60.30+- 2.5% 
59.21+- 3.0% 


1 aN 


11.48+- 8.9% 
41. 834=07.. 5% 
13.06+- 7.3% 
14.68+- 7.8% 
18.26+- 8.8% 
22.21+- 9.1% 
21.914-14.5% 
16.43+-17.9% 
6.52+-33.1% 
0.28+- 0.0% 
0.12+-16.7% 
0.12+- 8.3% 
0.12+- 8.3% 
0.12+-16.7% 
0.11+-18.2% 
0.12+-16.7% 
0.11+-18.2% 
0.12+-16.7% 
Q0.12+-16.7% 
0.12+-16.7% 
0.12+-16.7% 
0.12+-16.7% 
0.30+-83. 3% 
1.36+-58. 8% 
3.91+-29.4% 
6.381+-22.6% 
9.424+-22.1% 
12.27+-20.0% 
14.92+-18.1% 
16.77+-15.9% 
17.56+-13.8% 
17.21+-12. 5% 
16.61+-14.6% 
14.08+-13.7% 
9.39+- 9.9% 
8.414+- 9.2% 
7.80+-10.8% 
7.23+-10.8% 
6.85+-11.8% 
6.454+-13.6% 
5.71+-12.6% 
5.244+-11.8% 
4.76+-11.8% 
4.354+-12.0% 
4.05+-10.9% 
3.79+-10.3% 
3.52+- 9.9% 


L.TAN 
(AVG) 


0.150 
0.150 
0.163 
0.181 
0.225 
0.285 
0.333 
0.393 
0.351 
0.067 
0.039 
0.039 
0.039 
0.041 
0.041 
0.043 
0.042 
0.043 
0.042 
0.043 
0.043 
0.042 
0.068 
0.151 
0.253 
0.299 
0.305 
0.299 
0.287 
0.266 
0.243 
0.225 
0.212 
0.181 
0.125 
0.114 
0.108 
0.101 
0.097 
0.094 
0.087 
0.081 
0.075 
0.070 
0.066 
0.063 
0.059 


COND. 
MHO 


1.784 
1.841 
2.029 
2-281 
2.840 
3.456 


“3.405 


2-557 
1.012 
0.044 
0.019 
0.019 
0.019 
0.019 
0.017 
0.019 
0.017 
0.019 
0.019 
0.019 
0.019 
0.019 
0.047 
0.212 
0.608 
1.062 
1.464 
1.910 
Ze3A9 
2-608 
2.733 
2.677 
2.582 
2.193 
1.458 
1.307 
1.213 
1.125 
1.068 
1.005 
0.886 
0.817 
0.741 
0.679 
0.628 
0.589 
0.547 
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DIELECTRIC DATA FOR 1% DMSO IN WATER 


TEMP. 
DEG.C 


= NO 


pay fo) (eS) yey (op) (U0) 
6 
a ~) — o ~) ~) fO 


EP 


T5.02t- 1.8% 
T4.454+— 3.0% 
75.35+- 3.3% 
T4.26+- 5.0% 
60.56+- 6.9% 
43.53+-22.0% 
23<414-31.5% 
9.38+-29. 6% 
3.44+-16. 9% 


Si 02+- 8. 6% 
2.85+- 8.4% 
2-69t- 8.2% 


2374-11. 0% 
2- 544-10. 2% 
2.474+-12. 6% 
2.764+-12.7% 
2.95+-13.9% 
3.764+-15.7% 
5. 564-21. 4% 
8.99+-23.9% 
13.98+-24.7% 
20.79+-36. 5% 
29.72+-35. 3% 
41.15+-28. 1% 
53.03+-19.6% 
64.88+-17.4% 
T41. 574-13. 3% 
76.65+- 1.5% 
76.19+t- 2.0% 
75.73t+- 2.3% 
73.50+- 2.7% 
74.13t- 2.3% 
73.50+- 2.0% 
72.43+- 2.5% 
71.-57+- 2.1% 
67.33+- 5.3% 
65.39t- 3.4% 
64.32+- 2.8% 
63.50+- 2.8% 
62.044+- 2.9% 
61.29+- 3.1% 
60.19+- 2.9% 
59.11+- 3.0% 


EPP 


11,.24+- 7.7% 
12.81+-11.2% 
15.354+-13.3% 
18.22+-15.6% 
16.84+-13.8% 
13.80+-29.4% 
7.05+-50.2% 
1.60+-60. 6% 
0.33+-36.4% 
0.19+-26.3% 
0. 144+-21.4% 
Q0.12+-16.7% 
0.12+-25.0% 
0. 144+-21.4% 
0. 16+-18.8% 
0.22+-22.7% 
0.28+-21.4% 
0.38+-18.4% 
0.61+-27.9% 
1.204+-41.7% 
2o474+-45.7% 
4,63+-56.2% 
7.644+-51.6% 
10.80+-33.9% 
13.69+~30.3% 
16.13+-26.8% 
15.06+-35.6% 
14. 16+-16.0% 
13.08+-14.1% 
12.32+-11.7% 
10.02+-10.5% 
11.19+- 9.9% 
10.32+- 9.3% 
9.24+-10.0% 
8.50+- 8.0% 
7.39+-11.8% 
6.60+-12.7% 
6.084+-14.1% 
5.39+-11.5% 
5.07+-12.0% 
4.544+-12.1% 
4.214+-11.9% 
3.88+-11.9% 


L.TAN 
(AVG) 


0.150 
0.172 
0.204 
0.245 
0.278 
0.317 
0.301 
0.171 
0.096 
0.063 
0.049 
0.045 
0.051 
0.055 
0.065 
0.080 
0.0395 
0.101 
0.110 
0.133 
0.177 
0.223 
0.257 
0.262 
0.258 
0.249 
0.210 
0.185 
0.172 
0.163 
0.136 
On iat 
0. 140 
0.128 
0.119 
0.110 
0.101 
0.095 
0.085 
0.082 
0.074 
0.070 
0.066 


COND. 
MHO 


1.747 
1.982 
2.388 
2-834 
2.620 
2.149 
1.100 
0.249 
0.051 
0.030 
0.022 
0.019 
0.019 
0.022 
0.025 
0.034 
0.044 
0.059 
0.095 
0.187 
0.384 
0.723 
1.188 
1.678 
2na30 
2.208 
2-344 
2-205 
2.036 
1.916 
15568 
1.740 
1.602 
1.439 
Le ue) 
1.150 
1.024 
0.949 
0.836 
0.092 
0.704 
0.653 
0.603 
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DIELECTRIC DATA FOR 5% DMSO IN WATER 


TEMP. 
DEG.C 


EP 


7T4.904+- 2.0% 
73.12+- 2.0% 
T4.38+- 2.9% 
7T5.-4U8+- 3.2% 
72.06+- 5.9% 
54.27+-12.7% 
33.56+-22. 1% 
17.20+-29.3% 
7.724+-24.4% 
4.97+-34.0% 
3.78+=15. 1% 
3.67+-12.3% 
3.62+-13. 5% 
3.694t-14.1% 
4.104+-14., 1% 
4.50+-16.2% 
5.79+-16.4% 
7.97+-15. 3% 
10.63+-13. 0% 
15. 53+-20.9% 
20.94+-20. 1% 
27.72+-21. 1% 
36.394+-24.4% 
46.65+-21.8% 
59.16+-14.0% 
69.03+-12. 7% 
75.28t- 2.5% 
75.48+- 2.0% 
75.22+- 1.6% 
TW4.77+- 1.9% 
74.38t- 2.3% 
73.434+- 1.9% 
72. 74+- 2.4% 
71.944+- 2.3% 
71.08t- 2.3% 
67. O4+— ar 9% 
64.71+- 2.8% 
63-77+- 2.9% 
62.68+- 2.8% 
61.45+- 3.0% 
60.66+- 2.9% 
59.52+- 2.8% 
58.25+- 2.7% 


EPP 


11.98+-13.0% 
12.48+-10.3% 
14.78+-10.1% 
17.58+-13.4% 
21.05+-18.0% 
19.804+-27.4% 
14.25+-34.0% 
8.134+-34.4% 
3.244-44, 8% 
1.31+-71.8% 
0.58+-60.3% 
0.32+-43.7% 
0.424+-33.3% 
0.554-34.5% 
0.814+-33.3% 
1,.274+-34.6% 
2.104+-34.3% 
3.474+-31.7% 
5.25+-29.1% 
7.634+-37.1% 
9.91+-31.6% 
12.09+-30.8% 
14.434+-34.9% 
16.50+-31.1% 
18.104+-23.6% 
18.334+-24.6% 
17.444+-13.6% 
15.58+-11.9% 
14.00+-10.2% 
12.89+- 9.9% 
12.17+-11.7% 
10.77+- 9.3% 
9.98+-10.6% 
9.30+-10.1% 
8.844+- 9.7% 
7.874+-13.9% 
7.004+-10.4% 
6.45+-10. 5% 
5.88+-10.5% 
5.37+-10.6% 
5§.05+-10.5% 
4. 72+-10.4% 
4.35+-10.3% 


COND. 
MHO 


1.866 
1.942 
2-300 
2.733 
3.274 
5.070 
2.218 
1.263 
0.504 
0.204 
0.090 
0.050 
0.065 
0.085 
0.126 
0.197 
0.327 
0.540 
0.817 
1.188 
1.539 
1.879 
22243 
2.564 
2.815 
2.853 
2.714 
2-425 
2.180 
2.004 
1.891 
1.678 
1.552 
1.445 
1.376 
a2. 
1.087 
1.005 
Oe 
0.836 
0.785 
Wiel pete 
0.679 
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DIELECTRIC DATA FOR 10% DMSO IN WATER 


TEMP. 
DEG.C 


23.0 
15.5 


EP 


73.56+- 2.9% 
74.32+- 2.4% 
75.15+=- 2.4% 
75.67+—- 2.0% 
74.58+- 3.1% 
49.71+-21.3% 
26.474-27. 9% 
12.66+-46. 0% 
6.18+-20.7% 
5. 04+-14. 3% 
4.194+-12. 6% 
3.894-12. 3% 
4.03+-13.2% 
4.084+-13.7% 
4.53+-10.6% 
5.12+- 9.6% 
6.37+-17.1% 
8.484+-18. 8% 
12. 58+-35.3% 
22-054+-21. 8% 
33.20+-17.8% 
46.93+4+-14. 0% 
59.73+- 9.6% 
7O0.77+- 5.5% 
73.62+- 2.3% 
73.94+- 2.9% 
73.94+- 3.1% 
1a 18+- Ze 1% 
70.83+- 2.4% 
71.69+- 4.3% 
67-21t- 2.2% 
65.10+- 2.3% 
61.93+- 2.9% 
59.26+- 3.1% 


BEE 


11.994+-14.0% 
13.73+-14.2% 
16.004+-13.4% 
20.20+-14. 8% 
26.044+-15.9% 
23.30+-36.4% 
17.17+-39.4% 
§.38+-44.6% 
3.01+-29.9% 
1.354+-27.4% 
0.68+-26.5% 
0.42+-19.0% 
0.624+-19.4% 
0.834+-19.3% 
1.22+-18.9% 
1.884+-17.6% 
3.12+-22.1% 
5.104-21.0% 
8.304+-34. 5% 
13.87+-36.7% 
17.85+-36. 5% 
20.56+-30.8% 
21.134+-24.3% 
21.20+-19.1% 
18.994+-15.1% 
17.03+-15.2% 
14.49+-14.7% 
11.844+-13.0% 
10.664+-13.2% 


9.07+-11. 5% 
7.07+- 7.8% 
6.13+- 9.8% 
5.13+- 8.2% 
4.53+- 9.5% 


L. TAN 
(AVG) 


0.163 
0.185 
0.213 
0.267 
0.349 
0.469 
0.649 
0.662 
0.487 
0.268 
0.162 
0.108 
0.154 
0.203 
0.269 
0.367 
0.490 
0.601 
0.660 
0.629 
0.538 
0.438 
0.354 
0.300 
0.258 
0.230 
0.196 
0.164 
0.151 
Oa tee 
0.105 
0.094 
0.083 
0.076 


COND. 
MHO 


1.866 
2.136 
2.488 
3.142 
4.053 
3.625 
2-670 
1.301 
0.468 
0.210 
0.106 
0.065 
0.096 
0.129 
0.190 
0.292 
0.485 
0.792 
1.288 
2.155 
2.777 
3.198 
3.286 
3.299 
2.953 
2-652 
ZeZo0 
1.841 
1.659 
1.414 
1.100 
03955 
0.798 
0.704 
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DIELECTRIC DATA FOR 20% DMSO IN WATER 


TEMP. 
DEG.C 


23.0 
15.0 


EP 


7T4.83t+- 2.1% 
73-43t- 2.7% 
73.69+- 3.2% 
73.50+- 4.6% 
T2.-494+—- 4.2% 
67.50+- 8.0% 
30.43+4+-64.0% 

8.974+-21. 2% 

5.924+-14.9% 


5.-49+- 9.8% 
5.01+- 9.8% 
4.55+- 7.0% 
4.734+- 8.5% 
4.88+- 6.8% 
5.21+- 8.3% 


5.60+-10.5% 
6.95+-18.1% 
8.544+-17. 0% 
10.05+-10.6% 
10.57+- 7.3% 
53.90+-26.1% 
66.81+-15.3% 
72.99+- 7.4% 
75.48+- 3.2% 
75.09t—- 2.9% 
7T4.51t- 2.3% 
7T4.514+- 2.1% 
73.69+- 2.1% 
73.50+- 2.4% 
72.31t- 2.2% 
66.81t- 2.1% 
64.32+- 2.54 
61.13+- 3.0% 
58.65+- 3.4% 


EPP 


15.96+-10.3% 
17.67+-11.5% 
20.32+-13.1% 
26.06+-16.7% 
33.70+-15.9% 
4Q0.07+-19.4% 
23.93+-66. 0% 
8.074+-25. 8% 
3.90+-23.8% 
2.11+-20.9% 
1.194+-17.6% 
0.71+-18.3% 
1.014+-27.7% 
1.35+-30.4% 
1.944+-33.5% 
2.88+-36.5% 
4.58+-38.6% 
6.874+-33.2% 
8. 054-03 2% 
10.67+-19.3% 
30.58+-47.5% 
28.524+-29.3% 
26.19+-14.8% 
23.67+- 9.0% 
21.694+- 9.8% 
19.33+- 8.5% 
17.36+-11.0% 
15.01+-10.6% 
13.37+- 9.1% 
11.484+-10.5% 
8.86+-11.9% 
7.60+-13.2% 
6.24+-11.9% 
5.454+-11.7% 


L.TAN 
(AVG) 


0.213 
0.241 
0.276 
0.355 
0.465 
0.594 
0.786 
0.900 
0.659 
0.384 
0.238 
0.156 
0.214 
0.277 
0.372 
0.514 
0.659 
0.804 
0.839 
0.009 
0.567 
0.427 
0.359 
0.314 
0.289 
0.259 
0.233 
0.204 
0.182 
0.159 
0.133 
0.118 
0.102 
0.093 


COND. 
MHO 


2.482 
2.746 
3.160 
4.053 
52200 
6.233 
3.720 
PS257 
0.606 
0.328 
0.185 
0.111 
0.157. 
0.210 
0.302 
0.448 
0.710 
1.068 
1.313 
1.659 
4.756 
4.436 
4.072 
3.682 
3.374 
3.003 
26702 
2.337 
2.080 
1.784 
1.376 
1.787 
0.968 
0.848 
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DIELECTRIC DATA FOR 70% DMSO IN WATER 


TEMP. 
DEG.C 


22.0 
16.0 


EP 


51.45+- 5.4% 
31 .Dit=20. 2% 
26.53+-27.9% 
13.19+- 7.5% 
12.75+-12. 7% 
12.13+-49.7% 
7.99+-21. 7% 
7.-44+-10. 1% 
7.28+- 5.4% 


7.20+t- 5.6% 
6.78+- 3.1% 
6.83+- 5.9% 


7.20+- 4.9% 
7.244+-18.1% 
7.114+-22.2% 
9.63+-35.2% 
11.88+- 8.7% 
13.35+- 8.5% 
21.474+-25.9% 
24.83+-40.9% 
29.20+-43.7% 
33.03+-53.9% 
33.31+-55. 8% 
&6.77+-18.9% 
49OL4U14- 8.1% 
52.21+- 7.1% 
53.26+- 7.0% 
52-714+- 5.3% 
53.21+- 4.3% 
53.53+- 4.2% 
53.16+- 4.0% 
53.53+- 6.4% 
53.35+- 4.6% 
53.21+- 4.2% 
51.72+- 4.5% 
51.72+- 4.1% 
50.57+- 5.2% 
50.36+- 4.4% 


Eee 


26.414+-12.9% 
23.37+-34.1% 
19.04+-37.4% 
10.63+-16.9% 
9.05+-19.9% 
7.474+-44.2% 
W 372-90 4% 
3.05+-13.4% 
22.23+-11.7% 
1.69+-12.4% 


1.12+-19.8% 
2.104-12.9% 
2.904-99.0% 


3.904+-19.7% 

5.13+-20.9% 

7.68+-33.7% 
10.11+-17.5% 
12.03+-18.4% 
16.38+-42.2% 
17.74+-58.0% 
19.47+-56.3% 
20.72+-72.0% 
20.314+-75.7% 
24.79+-35.5% 
23.85+-24.2% 
23.18+-23.9% 
22.12+-23.3% 
22.99+-13.8% 
20.91+-12.5% 
19.32+-13.7% 
17.484+-13.4% 
15.86+-15.8% 
14.50+-13.0% 
13.14+-14.0% 
11.63+-13.7% 
10.69+-13.8% 

9.704+-15.4% 

8.95+-14.5% 


L.TAN 
(AVG) 


0.513 
0.623 
0.718 
0.806 
0.710 
0.616 
0.547 
0.410 
0.306 
0.235 
0.165 
0.307 
0.403 
0.539 
0.722 
0.798 
0.851 
0.901 
0.763 
0.714 
0.667 
0.627 
0.610 
0.530 
0.483 
0.444 
0.415 
0.436 
0.393 
0.361 
0.329 
0.296 
0.272 
0.247 
0.225 
0.207 
0.192 
0.178 


COND. 
MHO 


4.109 
3-638 
2.959 
1.602 
1.407 
1.162 
0.679 
0.474 
0.347 
0.263 
0.174 
0.327 
0.451 
0.606 
0.798 
1.194 
lao! 
1.872 
2.545 
2.758 
3.029 
3.223 
3.160 
3.858 
3.707 
3.607 
3.443 
3.975 
3.255 
3.003 
20421 
2.469 
20206 
2-042 
1.810 
1.665 
1.508 
1.395 
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DIELECTRIC DATA FOR 100% DMSO 


TEMP. 
DEG.C 


23.0 
Wa 
12.5 
14.0 
13.5 
10.0 
-4.0 


EP 


46.494+- 5.3% 
46.69+- 5.1% 
47.01+- 4.8% 
33.56+-10.8% 
20.794+-14. 1% 
9.18+-16.8% 
4.66+-13.5% 
3.46+-13. 0% 
3.32+-10. 5% 
3.18+-11.0% 
3.04+-11. 8% 
2-974-10.4% 
3.04+-11.3% 
3.064+-11. 1% 
2.-954+-11.2% 
3.25+-12.9% 
3.18+-14. 2% 
3.644+-13. 2% 
4.08+-12. 5% 
4.844+-16.1% 
5. 96+-16.4% 
7./24+-16. 3% 
8.46+-14.4% 
10. 12+-12. 3% 
12. 24+-13. 0% 
14.754+-14.2% 
18.20+-14. 5% 
21.694+-12.9% 
25.654+-14.0% 
30.79+-13.9% 
35.93+-13. 5% 
41.504+-11. 4% 
45.82+- 5.8% 
46.33+4- 5.0% 
W4,85+- 4.9% 
43.86+- 5.0% 
42.94+- 6.7% 
41.29+- 5.3% 
39.44+- 5.5% 
37.93+- 5.7% 


EPP 


13.06+-11.6% 
13.734+-11.4% 
15.35+-10.9% 
11.65+-22.6% 
6.38+-69.4% 
1.92+-38.5% 
0.60+-26.7% 
0.29+-27. 6% 
0.20+-20.0% 
0.16+-18.8% 
0.144+-14.3% 
0.13+-15.4% 
0. 15+-20.0% 
0.16+-18.8% 
0.26+- 0.0% 
0.21+-19.0% 
Q0.25+-24.0% 
0.344+-20.6% 
0. 46+-19.6% 
0.65+-26.2% 
0.93+-28.0% 
1.29+-30.2% 
1.75+-28.0% 
20494-24.5% 
3.53+-25.2% 
4.83+-22.2% 
6.414+-21.2% 
7.76+-18.7% 
9.42+-21.5% 
11.15+-20.9% 
12.53+-19.6% 
13.62+-21.1% 
13.63+-16.2% 
119). 324 11, 
10.89+-13.9% 
9.484+-14.9% 
8.00+-17.3% 
6.79+-15.9% 
5.72+-12.1% 
5.03+-13.7% 


L.TAN 
(AVG) 


0.281 
0.294 
0.327 
0.347 
0. 307 
0.209 
0.129 
0.084 
0.060 
0.050 
0.046 
0.044 
0.049 
0.052 
0.088 
0.065 
0.079 
0.093 
0.1413 
0.134 
0.156 
0.178 
0.207 
0.246 
0.288 
0.327 
0.352 
0.358 
0.367 
0.362 
0.349 
0.328 
0.297 
0.266 
0.243 
0.216 
0.186 
0.164 
0.145 
0.133 


COND. 
MHO 


2.029 
2.136 
2.388 
1.810 
0.993 
0.298 
0.094 
0.045 
0.031 
0.025 
0.022 
0.020 
0.023 


20025 


0.040 
0.033 
0.039 
0.053 
0.072 
0.101 
0.145 
0.200 
0.272 
0.387 
0.549 
0.754 
0.999 
1.206 
1.464 
1.734 
1.948 
2-117 
2.117 
1.916 
1.696 
1.477 
1.244 
1.056 
0.892 
0.785 


237 


i 


2588 
perereregit se poreeee eras 


rv.2@ —*t2 TE 


o 
: 


DIELECTRIC DATA FOR 5% DEXTRAN IN WATER 


TEMP. 
DEG.C 


ee 
Vis 
11.3 


EP 


T4.324+- 2.9% 
75.15+- 5.2% 
T5.414+- 5.7% 
T4.454+- 5.0% 
69.56+-14. 2% 
43.574+-14. 3% 
24.964-17. 4% 
11.244+-20.4% 
3.67+-10. 6% 


3.06+- 8.8% 
3.02+- 8.6% 
2.97+- 6.7% 
2.76t+- 7.2% 
2.71+-10. 3% 


2.87+-10. 8% 
2.78+-11. 2% 
2-90+-10.3% 
S8O0U+—-112 2% 
34094-12235 
3.62+-11. 3% 
5.43+-15.7% 
8.914+-15.2% 
14.45+-17.2% 
21.30+-18. 5% 
31.89+-21. 1% 
41.22+-13.1% 
54.65+-11.9% 
64.60t+- 8.3% 
75.67+- 9.1% 
77.31+- 4.6% 
76.194+- 4.2% 
75.80+- 4.1% 
7T5-48+- 4.0% 
75.41t- 4.1% 


EPP 


14.134+-14.9% 
Voce etm oo eS 
18.08+-17.9% 
21.214+-14.9% 
24.304+-30.1% 
16.82+-20.9% 
9.79+-28.8% 
22.214+-72.4% 
0.25+-36.0% 
0. 144+-14.3% 
0.134+-15.4% 
0.12+-16.7% 
0.12+-16.7% 
0.12+-16.7% 
0.134+-15.4% 
0.144-21.4% 
0.15+-20.0% 
0. 164-12.5% 
0.18+-16.7% 
0.33+-18.2% 
0.714+-15.5% 
1.424+-23.9% 
3.32+-34.6% 
6.18+-36.1% 
10.40+-40.2% 
13.25+-24.9% 
16.32+-17.6% 
18.73+-16.6% 
21.23+-15.1% 
19.86+-13.5% 
17.82+-11.7% 
16.67+-11.1% 
15.20+- 9.7% 
14.09+-12.6% 


L.TAN 
(AVG) 


0.190 
0.203 
0.240 
0.285 
0.349 
0.386 
0.392 
0.197 
0.068 
0.046 
0.043 
0.9040 
0.043 
0.044 
0.045 
0.050 
0.052 
03053 
0.058 
0.091 
0.131 
0.159 
0.2390 
0.290 
0.326 
ees 81 
0.299 
0.290 
0.281 
0.257 
0.234 
Us22u 
0.201 
0.187 


COND. 
MHO 


2.199 
2.369 
2.815 
a fey fe) 
3.782 
2.614 
1.521 
0.344 
0.039 
0.022 
0.020 
0.019 
0.619 
0.019 
0.020 
0.022 
0.023 
0.025 
0.028 
0.051 
0.111 
0.221 
0.516 
0.961 
1.615 
2.061 
2.538 
Pegs Ht fs 
3.305 
3.091 
Zot 
rey abs) 
2-362 
Pa REE 
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DIELECTRIC DATA FOR 0.1% NACL IN 10% DMSO 


TEMP. 
DEG.C 


EP 


73.75+- 5.7% 
76.19+- 4.1% 
76.78+- Yiug 
76.98+—- 3.5% 
68.144+-17. 2% 
44.514+-23.0% 
22.-82+-24. 0% 
10.104+-22. 4% 
6.494+-13.3% 


4.30+-17.2% 
3.644+- 9.9% 
3.464+-10. 4% 
3.13+- 9.9% 
3.11+- 9.6% 
3.50+- 9.7% 
3.644+-10.2% 
4.41+- 9.1% 


5.43+-13.4% 
7.72+-18. 8% 
11.72+-25.7% 
17.294+-35. 5% 
28.644+-31.9% 
44.014+-21.4% 
59.06+-17.9% 
69.87+-10.8% 
74.64+— 7.1% 
74.90+- 3.5% 
T4.324- 2.4% 
74.58+- 3.1% 
73.754+- 2.2% 
72.99+- 1.8% 


ERE 


14.324+-11.5% 
15.79+-10.3% 
18.41+-11. 0% 
21.77+-13. 6% 
23.91+-31.5% 
16. 24+-39.2% 
10.97+-30. 5% 
4.59+-26.1% 
1.974+-21.8% 
G.81+-24.7% 
O.404-20.0% 
0.29+-20.7% 
0.23+-21.7% 
0.32+-25.0% 
0.444-27.3% 
0.62+-27.4% 
0.99+-29.3% 
1.63+-33.7% 
2.90+-37.9% 
5.01+-39.7% 
7.87+-45.7% 
12.91+-46.5% 
16.91+-38.7% 
19.40+-34.5% 
19.35+-29.6% 
18.63+-14.9% 
15.88+-18.0% 
14.544+-17.3% 
13.18+-18.7% 
12.02+-17.2% 
11.15+-15.7% 


L.TAN 
(AVG) 


0.194 
0.207 
0.240 
0.283 
0.351 
0.432 
0.481 
0.454 
0.304 
0.188 
0.110 
0.084 
0.073 
0.103 
0.126 
0.170 
0.224 
0.300 
0.376 
0.427 
0.455 
0.451 
0.384 
0.328 
On 207 
0.250 
0.212 
0.196 
Ors ay, 
0.163 
0.153 


COND. 
MHO 


2.224 
2-457 
2.865 
3.387 
3.720 
2-991 
1.709 
0.716 
0.307 
0.126 
0.062 
0.045 
0.036 
0.050 
0.068 
0.096 
0.154 
0.253 
0.451 
0.779 
1.225 
2-011 
2.633 
3.016 
3.010 
2-897 
2.469 
2-262 
2-048 
1.872 
1.734 
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DIELECTRIC DATA FOR 1.0% NACL IN 10% DMSO 


Lh ays Wc 
DEG.C 


EP 


73.62+- 2.3% 
73.56+- 2.9% 
73.62+- 3.2% 
72.06+- 4.7% 
62.-414+- G.4% 
40.83+-14.9% 
18.62+-20. 2% 
10.574+-14.1% 
7364 1026% 


4.66+- 5.2% 
3.85+- 8.6% 
3.76+- 6.1% 
3.444+-— 6.1% 
3.444+- 8.1% 
3.87t- 7.5% 
G.21+- 5.0% 
5.17t- 8.1% 
7.16+-11. 2% 
9.99+-14.7% 


12.28+-13.4% 
18.65+-19.0% 
28. 864+-18.4% 
41.78+-15. 8% 
60.71+- 9.3% 
69.81+- 6.8% 
THU.83+- 3.0% 
T4.514+- 2.7% 
73.56+- 3.1% 
73.12+- 2.1% 
72.81+- 2.3% 
72.68+- 2.3% 


EPP 


21.39+- 9.8% 
21.004+- 7.3% 
22.57+- 7.4% 
25.074+-— 1028505 
25.994+-14.6% 
20.244+-23.8% 
10.58+-23.8% 
52394-21208 
2.-65+-23.4% 
1.144+-15.8% 
0.56+-14.3% 
0.364+-11.1% 
0.30+-10.0% 
0.42+-14.3% 
0.61+-14.8% 
0.98+-15.3% 
1.66+-20.5% 
3.27+-42.5% 
4.65+-18.5% 
6.484-17.7% 
10.28+-20.9% 
14.92+-25.1% 
19.42+-23.6% 
25.07+-20.9% 
25.934+-15.1% 
24.72+-10.1% 
23.01+- 8.9% 
21.70+- 8.6% 
21.08+- 6.5% 
20.67+- 6.3% 
20.33t- 7.2% 


L.TAN 
(AVG) 


0.291 
0.285 
0.307 
0.348 
0.416 
0.496 
0.568 
0.510 
0.360 
0.245 
0.1745 
0.096 
0.087 
6.122 
0.158 
0.233 
0.321 
0.457 
0.465 
0.528 
0.551 
0.517 
0.465 
0.413 
0:°371 
0.330 
0.309 
0.295 
0.288 
0.284 
Ono 


COND. 
MHO 


3.324 
3.267 
3.512 
B59 02 
4.040 
3.148 
1.646 
0.836 
0.412 
0.177 
0.087 
0.056 
0.047 
0.065 
0.095 
0.153 
0.258 
0.508 
0.723 
1.005 
1.596 
7A HB) 
sree 
32902 
4.034 
3.845 
3.581 
3.374 
3.280 
3.217 
3.242 
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DIELECTRIC DATA FOR 5% DEXTRAN IN 10% DMSO 


TEMP. 
DEG.C 


22.1 
20.6 
15. 9 
8.3 
30-2 
=i. 7 


EP 


67.-444+- 2.8% 
71.87+- 3.0% 
72.06+- 3.2% 
73.37+- 4.3% 
73.37+- 4.5% 
68. 26+-12. 6% 
39.20+-36.4% 
266 324-35). 1% 
11.53+-44. 3% 
7.07+-36.9% 
4.59+-30.3% 
4.034+-14.6% 
3.55+-10. 1% 
3.46+- 9.0% 
3.78+- 7.4% 
G.O1+- 7.7% 
Wo77+- 7.1% 
5.90+-11. 5% 
8.44+-26. 3% 
11.81+-22.1% 
16.64+-24.9% 
25.91+-26.9% 
39.64+-38. 8% 
50.10+-22.4% 
63.28+-14.4% 
67.444+- 7.1% 
68.38+t- 5.7% 
68.26+- 5.4% 
68.20+- 4.5% 
68.32+- 5.0% 
68.03+- 4.7% 
67.79+- 4.4% 


EPP 


13.064+-12.5% 
13.66+-12.7% 
14.73+-12.9% 
17. 46+-15. 6% 
21.39+-19. 6% 
22.61+-35.8% 
16.38+-75.3% 
11.74+-46.7% 
5.09+-77.6% 
2.07+-76.8% 
0.804+-57.5% 
0.41+-36.6% 
0.25+-24.0% 
0.27+-22.2% 
0.46+-13.0% 
0.66+-15.2% 
1.07+-15.9% 
1.83+-18.0% 
3.34+-30.8% 
52514+-25.0% 
8.744+-25.5% 
12.93+-34.4% 
17.35+-52.1% 
19.55+-38.6% 
21.534+-29.7% 
19.84+-17.4% 
18.09+-14.3% 
16.60+-11.8% 
15.454+- 9.3% 
13.71+-10.2% 
12.-244-12.4% 
11.364+-14.9% 


L.TAN 
(AVG) 


0.194 
0.190 
0.204 
0.238 
0.292 
0.331 
0.418 
0.503 
0.441 
0.293 
0.174 
0.102 
0.070 
0.078 
Oe 122 
0.165 
0.224 
0.310 
0.396 
0.467 
0.525 
0.499 
0.438 
0.390 
0.340 
0.294 
205 
0.243 
0.227 
0.201 
0.180 
0.168 


COND. 
MHO 


2.029 
2.124 
ae 293 
2-714 
3.324 
= Pape va) 
2.545 
1.828 
0.792 
0.322 
0.124 
0.063 
0.039 
0.042 
0.072 
0.102 
0.167 
0.285 
0.520 
0.855 
1.357 
2.011 
2-695 
3.041 
3.349 
3.085 
2-815 
Zeooe 
2-400 
2-130 
1.904 
1.766 
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DIELECTRIC DATA FOR DMSO IN WATER AT 5 C 


CONC. 


EP 


80.00+- 2.8% 
75.00+- 3.6% 
73.444+- 3.3% 
71.384- 3.3% 
62.47+- 2.9% 
50. 18+-10.6% 
38.20+-13.5% 
27.00+=- 8.2% 
32.50+- 9.3% 
43.34+- 5.6% 
8.42+- 1.8% 
8.424+- 1.8% 


BEEP 


14.95+- 8.4% 
18.77+- 9.7% 
23.65+- 7.6% 
30.03+-10.1% 
34.09+- 8.4% 
32.944+-17.2% 
28.00+-19.1% 
20.844+-13.9% 
24.50+-16.7% 
25.944+-12.6% 
25.94+-12.6% 

12 81t— 92225 


DIELECTRIC DATA FOR DMSO IN WATER AT 


CONC. 


EP 


76.58t- 1.9% 
73.24+- 2.0% 
72.06t+- 3.5% 
70.53+- 2.6% 
67.21+- 3.1% 
65.22t- 2.8% 
61.03+- 3.9% 
57.95+- 5.8% 
55.60+- 4.2% 
52.80+- 4.2% 
4u7.50+- 5.7% 


EPP 


11.48+- 8.9% 
11.99+- 8.4% 
13.91+t- 9.8% 
17.31+- 8.8% 
20.744+- 9.6% 
24.97+-11.8% 
27.63+-13.2% 
28.45+-14.5% 
25.93+-11.1% 
20.40+- 9.8% 
13.42+-10.1% 


L.TAN 
(AVG) 


0.187 
0.250 
0.328 
0.322 
0.546 
0.656 
0.733 
0.772 
0.754 
0.598 
0.599 
0.215 


COND. 
MHO 


Zetno 
Auaee 
3.676 
4.668 
5-303 
et 
4.355 
3.239 
3.813 
4.034 
0.281 
0.281 


COND. 
MHO 


1.762 
1.862 
2.161 
2-689 
s2223 
Joos 
4.298 
4.423 
4.034 
Sods 
2.086 
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DIELECTRIC DATA FOR DMSO IN WATER AT 50 C 


CONC. 


EP 


66. 64+- 
66.35+- 
65.89+- 
65295+— 
64.55+- 
61.40+- 
55.03+- 
54.23+- 
52.44+- 
50. 14+- 
46.77+- 


3.6% 
3.3% 
3.4% 
3.3% 
3.3% 
5.3% 
5.4% 
7.0% 
52 2% 
4.8% 
6. 5% 


EPP 


7.58+-12.1% 
8.79+-12.1% 
11.12+- 9.8% 
13.044+-11.5% 
16.12+-17.9% 
19.41+-13.3% 
18.97+-11.1% 
18.35+-20.3% 
17.594+-12.4% 
14.404+-12.2% 
10.70+-22.1% 


L.TAN 
(AVG) 


0.114 
0.132 
0.169 
0.199 
0.250 
0.316 
0.345 
0.338 
0.335 
0.287 
0.229 


COND. 
MHO 


1.181 
1.370 
1.728 
2.029 
Paps) OH | 
3.016 
2.953 
2.853 
Za 
2.237 
1.665 
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